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NOIILINEAR OSCILLATIONS OF HAPMONIGALLY EXCITED 
AUTOPARAMBTRIC STSTMS 
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by 
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As the title suggests, this thesis presents a 
study on nonlinear oscillations of harmonically excited 
autqparametric systems* Two different types of systems 
are investigated* The first one has two degrees of 
freedom* It comprises of a spring - mass - dashpot pri- 
mary system and a secondary system which is a damped 
rigid pendultun hinged to the primary mass* The second 
type of systan, having three degrees of freedom, is very 
similar to the first one except that in this case the 
secondary consists of an elastic pendulum. In both cases 
the primary mass is excited by a harmonic force of cons- 
tant amplitude* The linear natural frequencies of both 
the ^sterns are adjusted so as to cause internal resonance. 
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The work is hroadly divided in five sections* 

The first four sections are devoted to the study of sys- 
tems of the first type. The last section deals ^?ith only 
the absorber action of the secondary in both types of 
systems (i.e.^with rigid and elastic pendulm)* 

The first part deals with the harmonic responses 
under small excitation. Two types of restoring couples on 
the pendulum are considered. In the first case the resto- 
ring couJ)le is entirely due to gravity field, when the 
pendulum can execute continuous revolutions. In the second 
case a linear torsional spring, preventing continuous 
revolutions of the pendulim, provides the entire restoring 
couple. VJith a harmonic forcing function, the steady state 
responses are also expected to be harmonic. As such the 
method of harmonic balance is used. Eie stability charac- 
teristics are studied by perturbing the approximate solu- 
tions. >Jhenever necessary , the analytical results are 
successively refined with the help of higher order approxi- 
mations by including higher order nonlinearities. The 
approximate results are verified by numerically integrating 
the equations of motion. 

With both types of pendulums (spring ccntrolled 
and gravity controlled) it is observed that for an excita- 
tion amplitude greater than a critical value (which 
depends on other system parameters) the responses no 
longer remain harmonic in certain frequency ranges* 



Integrated resiiLts show thi.t v?ith the spring controlled 
pendiiLuin the responses still remain periodic, whereas, 
the gravity controlled pendiiLum may result in nonperio- 
dic responses* For the system having a spring controlled 
pendulum, the stable harmonic motion bifurcates to one 
unstable harmonic motion and another stable amplitude 
modulated periodic state. This bifurcation is found to 
be associated with a jump in the amplitude values. Eest 
of the thesis is devoted only to the system with a gravity 
controlled pendulum. 

Higher values of the excitation parameter are 
dealt with in the second part of the thesis, in the re- 
gions of unstable harmonic solutions, numerical integra- 
tions are carried out to study the nature of the responses. 
Depending on the frequency, the following two types bran- 
ching of the stable harmonic motion are obtained: 

i) With a moderate excitation the stable harmonic 
motion branches out to one (or more) unstable 
harmonic state (s) and one (or more) amplitude 
modulated periodic motion(s). The amplitude 
modulated states (depending on the initial con- 
ditions) are weakly stable and the modulation 
period can vary over a wide range. In this bran- 
ching the amplitudes do not increase sharply. 
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ii) The stable hamionic motion branches out to one 

(or more) unstable harmonic state (s) and another 
chaotic state. This branching is associated 
with sharp increase in the amplitudes. 

The approximate method is not applicable in the 
chaotic regime and the results are obtained by only nume- 
rical Integration, which do not converge and the responses 
are completely nonperiodic. Thus the responses, though 
remaining bounded, can not be described uniquely. 

In the tliird part of this work statistical des- 
cription of the nonperiodic responses is attempted with 
the proposition that in such cases the actual physical 
experiment would yield randan responses and the integra- 
ted results represent these. Mean square values and PSD’s 
are obtained by using different step sizes of integration 
and initial conditions. These values are found not to 
differ much and as such are considered as meaningful des- 
criptions of the responses. The nature of the PSD func- 
tions for the primary and the secondaiy is investigated^ 
no attempt, hov/ever, is made to correlate the statistical 
descriptions with the system parameter combinations. 

In the fourth part of tlie thesis sane experimen- 
tal results are reported. The major objective of these 
experiments is to verify qualitatively the analytical 
results presented earlier. In the experiment, a harmonic 
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base excitation, instead of a force excitation, is given 
to the primaiy system* This does not alter the basic 
characteristics of the system investigated earlier. It 
is seen that for certain conbinations of parameters the 
system responses indeed become random and for these para- 
meter combinations the numerical integrations also do 
not converge. 

The vibration neutralizing effect of the pendu- 
lum is studied in the last part. The rigid pendulum in 
iiie two degree- of-freedcm system is shovm to exhibit 
limited absorber action. The range and scope of this 
absorber are discussed. The optimum value of the absor- 
ber damping is found to depend upon the level of excita- 
tion amplitude. lastly, the absorber action of a damped 
elastic pendulum is studied briefly • It is shown that 
with a nonop timal auxiliary mass damper, the peaks in the 
responses may be suppressed by providing the secondary 
with an additional degree of freedom in rotation. 



CHAPTER I 


B^TRODUCTION 

1.1 INTRODUCTION AND LITERATURE SURVEY 

Vibratory systems, in practice, are all nonlinear, 
but within certain limits they can, in many ways, be depic- 
ted in approximate fashion by a linear autonomous system 
for which the complete solution can be obtained in closed 
form. Unlike linear systems, nonlinear systems and non- 
autonanous linear systems do not lend themselves to general 
solutions. Moreover, many intricate features of the 
nonlinear systems are lost in the process of linearization. 
Systems of special interest are those possessing periodic 
solutions. Since the qualitative and the quantitative 
characteristics of such systems can not be obtained in 
totality, some approximate methods are called for to analyse 
these. Because of lack of generalized analysis, different 
types of nonlinear equations exhibit different interesting 
features. A special class of nonlinear systems includes 
parametric and autoparametric systems. 

Parametric vibration refers to the oscillatory 
motion that occurs in a system as a result of time-dependent 
(usually periodic) coefficients in the governing differen- 
tial equation Cl.e. nonautonomous equations). In such 
cases the time dependence is explicit, which implies an 
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external energy source (pump). Parametric oscillations 
may stabilize a system v/hich. is statically unstable or may 
destabilize a statically stable system. The standard second 
order hcciogeneoiis linear equations "with periodically time 
varying coefficients have the forra of Mathieu - Hill 
equations [1, 2], According to linear theory, the oscilla- 
tions in the regions of instability increase unboundedly 
and rapidly • These growths are checked due to nonlineari- 
ties of the system and the oscillations may beccme periodic. 
Various types of nonlinearities in a dynamical system are 
listed in reference [ 3] . 

For a two degree- of-freedcm system, the auto- 
parametric oscillations arise due to nonlinear coupling 
between the governing equations rather than the explicit time- 
dependence of the coefficients. This type of system is 
termed an autoparametric system (because no external pump- 
ing is required) or a system with internal resonance 
(because certain combinations of linear natural frequencies 
of the two modes cause resonance). The term autoparametric 
is more appropriate because unlike other resonances (which 
are observed at discrete frequencies only) these systems 
exhibit resonance over a frequency range - a characteris- 
tic of the parametric oscillations, 

A system, besides having parametric (or auto- 
parametric) oscillations, may be subjected to external 
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inputs appearing as nonhcmoceneous terms in the equations. 

An excellent review of parametric vibrations is given in 
reference - s] . Other useful references may be found 
in [9, 10] * 

In the present \fork the autoparametrlc vibrations 
of a simple and an elastic pendulum are studied. A simple 
pendulun \;ith an oscillating support has been used exten- 
sively to explain parametric instability. The response 
characteristics of this system have also been studied in 
detail, to quote a fei; references, Skalak and Yarymovytch 
[11], Caughey [12, 13], and Struble [*14, 1^] studied the 
steaifiy state responses of a pendulum v/ith a vertically or 
horizontally oscillating support. In these analyses only 
the first nonlinear term was retained. Chester [l6] inves- 
tiga,ted the same system with both horizontal and vertical 
motions of the support. 

Hie pendulum as a part of a two-degree- of-freedom 
autoparametrlc system has also been studied. The autopara- 
metric instability and energy transfer associated with 
coupled differential equations of a beam-pendulum system 
have been investigated [l7 - 19j • A similar system with a 
small external excitation \-ja.s studied in [20] . Free osci- 
llations of an elastic pendulum, with the linear natural 
frequencies in the longitudinal and angular motions adjusted 
so as to cause internal resonance, were investigated in 



[ 21 - 2^], Ryland and Meirovitch[ 25] determined the 
stability boundaries of a swinging spring with an oscilla- 
ting support. 

A part of the system, under investigation in this 
work, is a hinged pendulum capable of executing large 
motions. With first order approximation, the equations of 
motion have quadratic nonlinearities. Such equations have 
been studied in detail under the conditions of internal 
resonance. It is well known that with harmonic excitation 
of these systems, the stable periodic solutions may branch 
out to one unstable and two stable periodic solutions; this 
is associated with jumps in the amplitude values. It was 
shown in[ 26, 9] that if the excitation frequency is in the 
neighbourhood of the lower natural frequency (which is half 
of the other), another branching of the stable periodic 
solution may take place without any jump. In this case one 
unstable periodic solution and another stable amplitude 
modulated motion are obtained. Such phenomena have also 
been observed by Miles [27] in a system with two equal 
linear natural frequencies and with cubic nonlinearities, 
and by Yamamoto and Yasuda [ 28] in a system with one linear 
natural frequency twice the other and having quadratic and 
cubic nonlinearities. Sethna and Bajaj [29] derived the 
conditions on the system parameters at which the stable 
periodic motion bifurcates into an unstable periodic motion 
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and an amplitude modulated notion at the excitation fre- 
quency. In the course of present study also, such 
amplitude modulated motions are found but vith the differ- 
ence that there nay exist more than one set of amplitudes 
at which the modulation occurs and that these states are 
weakly stable. Moreover, another branching of the stable 
periodic states is observed where besides one (or more) 
unstable periodic state (s) a nonperiodic state is obtained. 
The latter state reveals almost randm motion of the system. 

The deterministic systems , with simple mathema- 
tical models, producing random or chaotic motions have ’long 
been observed in various biological, economic, social and 
physical sciences and have been reviewed by May [ 30] • 

Lorenz [ 31 ] called this general phenomenon the 'butterfly 
effect" : even if the atmosphere could be described by a 
deterministic model, the fluttering of a butterfly ^s wings 
could alter the initial conditions, and thus (in the chaotic 
regime) alter the long term predictions. Holmes [32] 
studied a nonlinear osci3-lator- for which complicated non- 
periodic motions arise. He called these motions as strange 
attractors, studied their structure by means of the 
Poincare map and developed a necessary criterion for 
strange attractor motions lAhich determines the minimum 

+ References [32 - 3^1 were noticed towards the end of 
this thesis* 
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forcing amplitude as a function of forcing frequency. As a 
mechanical example, Moon and Holmes [ 33 j examined the non- 
linear forced vibrations of a cantilever bean which is 
buckled by magnetic forces. Moon [ 3^] later presented an 
experimental criterion on the forcing amplitude and frequency 
for strange attractor motions and developed a heuristic 
formula. The basic system in [ 32 - 3^] had two stable and 
one unstable static equilibrium positions. ?or fixed damp- 
ing and frequency, the motion for small forcing amplitudes 
was periodic (about either of the stable fixed point) but 
■became chaotic for larger forcing amplitudes, with the 
solution jumping from one equilibrium to the other. For 
sufficiently large forces the oscillations encompassed all 
three fixed points of the unforced system resulting in 
periodic motions. The present work shows a possibility of 
such chaotic motions with a tvjo degree-of -freedom system 
where one mode is parametrically excited and the unforced 
system has only one static equilibrium position. 

Lastly a few applications of the parametric pheno- 
menon have also' been reported in the literature. In 
electrical devices this phenomenon has been exploited con- 
siderably [ 35 ] . Rodgers [36, 37] suggested the mechanical 
amplifiers based on the parametric phenomenon, i;hich was 
further studied by Eisinger and Merchant [38, 39] • Haxton 
and Barr [ 20 J suggested an autoparametric vibration absorber 
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(^/liich is rather a neutralizer). Curiousl;^ the practica- 
bility of such absorbers \‘JSls not explored further. Unlike 
t]ie conventional auxiliary mass damper which gives linear 
response characteristics, the auxiliary system (which is a 
cantilever beam with a concentrated mass at the free end) 
considered in |20] was excited because of nonlinear coupling. 
The range of frequency over which the cantilever beam osci- 
llates -wa-s determined from the instability zones of its 
static equilibrium position. The response characteristics 
03 ? such systems obviously depend on the force amplitude. 

In the present work a similar system (with a hinged pendulum 
instead of a cantilever beam) is studied in detail for 
different levels of excitations (harmonic) and is sliovni to 
exhibit many interesting features. 

Two ma3 or approaches to the solution of nonlinear 

differential equations are (i) the qualitative methods 
and (ii) the quantitative methods. Only the latter approach 

is used since we are Interested in the quantitative analysis 
as well. Some of the popular methods are methods of avera- 
ging ■ duo to Krylov - Bogoliubov - Mitropolski [^oj and 
Struble [4-l] , method of miTLtiple scales [42] and method 
of harmonic balance [3] • These methods have found inte- 
resting and wide applications to study various phenomena 
in nonlinear differential equations with or w’ithout 
excitation. The method of harmonic balance is not very 
suitable for the study of free oscillations where a priori 
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Imowledge of the solution is not availahle. For forced 
oscillations with harnonic inputs, hoiro^'or, this nethod is 
more straightforward to study the nature and amplitudes of 
the harmonic responses* Most of the results in this thesis 
are obtained using the method of harmonic balance. 

1*2 OBJECTIVE AND SCOPE OF THE PEESMT WOEK 

Nonlinear systems are mostly analysed for small 
oscillations. The objective of the present work is to 
study an autoparametric system under harmonic excitation 
with no restrictions on the amplitudes. The systtan comprises 
of a primary, which is a spring - mass - dashpot system, and 
a secondary, which is a damped rigid pendulum hinged to the 
primary mass. The characteristics of the spring and the 
dashpot s are assumed to be linear* The primaiy mass is 
excited by a harmonic force and its linear natural frequency 
is taken as twice that of the pendulum so as to cause inter- 
nal resonance* 

In Chapter 2 the existence and characteristics 
of the harmonic responses of this system are studied. 
Approximate results are obtained by using the method of 
harmonic balance. The accuracy of the solutions are checked 
by considering different orders of approximations. Stabi- 
lity of the haimonic motion is checked by perturolng the 

To check the validity of the approximate results 


solutions • 



9 


the equations of motion are numerically integrated* Further 
to stuc3y the effect of allowing the pendulum to make com- 
plete revolutions, two types of restoring couples on the 
pendulum are considered* In the first case the restoring 
couple is entirely due to gravity field, when the pendulum 
can execute continuous revolutions* In the second case a 
linear torsional spring, preventing continuous revolutions 
of the pendulum, provides the entire restoring couple* In 
both cases it is observed that for an excitation amplitude 
greater than a critical value (idiich depends on other system 
parameters) the responses no longer remain harmonic in 
certain frequency ranges. From the integrated results it 
is seen that the responses with the spring controlled 
pendulum still remain periodic, whereas, the gravity con- 
trolled pendulum results in nonperiodic responses* 

For excitation amplitudes greater than this 
critical value, further analysis is carried out in Chapter 3 
with the gravity controlled pendulim only. The major 
objective in this chapter is to determine the regions of 
stable harmonic solutions and the b^aviour of the respcmsos 
when the harmonic scCLuticns are unstable. To this end, 
the regions of stable and unstable harmonic solutions are 
obtained by the approximate method, ty plotting these 
regions in the forcing amplitude and frequency plane, it is 
observed that there are two separate regions where the 
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harmonic solutions are unstahlo. Numerical integrations 
reveal that only in one of these regions, and for a limited 
range of forcing amplitude, frequency and initial conditions, 
the response may be amplitude modulated. For other cases 
(where harmonic solutions are unstable) the responses are 
seen to be completely nonperiodic. For this reason the 
approximate solutions are not attempted in these ranges of 
unstable solutions. Father the empiiasis is placed on the 
results obtained by numerical integrations. 

For these cases, the integrated responses are seen 
to be nonperiodic and nonconvergent . Based on physical 
arguments, an attempt is made to Interpret tliese nonconver- 
gent results. It is proposed that in such cases the actual 
physical experiment would yield random responses and that 
the integrated results represent these. With this consi- 
deration, the statistical analysis of the integrated respon- 
ses is undertaken in Chapter 4-. We try to establish the 
randomness of the responses tiircwigh autocorrelation and 
power spectral density functions. For obtaining these 
statistical quantities the hypothesis of stationarity is 
assumed. An attempt is also made to justify this hypothesis 
by performing run tests. 

The statistical analysis in Chapter 4 is carried 
out with the expectation that the actual physical system 
(for certain combinations of parameters) would also result 



11 


in nonr op reducible and nonperiodic (i.e. random) responses* 
With this point of view an experiment is performed and the 
results are presented in Chapter 5* ob;iective of the 

experiment is only to verify the existence of random res- 
ponses of the deterministic system with a harmonic excita- 
tion, No statistical analysis of the experimental results 
is attempted. 

In Chapter 6 the objective is to stuc3y the vibra- 
tion - neutralizing effect of the pendulum. The pendulum, 
as a parametrically excited seccxidary system, is indeed 
found to neutralize the vibrations of the primaiy mass, which 
is acted upon by an external force. The range and scope of 
this absorber action arc discussed. Further, based on these 
results, the absorber action of a damped elastic pendulum 
is studied briefly. 



CHAPTER 2 


RESPONSE OF A TWO DEGREE - OF - FRiEDCM AUTOPARAMETRIC 
SISTEM TO HARMONIC EXCITATION 

2.1 INTRODUCTION 

In this chapter the forced oscillation of a two 
degree - of - freedon autoparametric systeu is studied. 

The system under investigation is shown in Fig. 2.1. The 
excitation is considered to ho hamonic and the range of the 
amplitude is limited so as to produce a harmonic response 
of the system. The system comprises of a primary, consist- 
ing of a linear spring - mass - dashpot and a secondary, 
consisting of a damped rigid pendulum. The primary mass is 
excited hy a harmonic force. The linear natural frequencies 
of the primary and the secondary are adjusted so as to cause 
internal resonance. The range of excitation frequency 
considered is such that the pendulum is parametrically 
excited. Two types of restoring couples on the pendulum 
are Investigated. In the first case, the restoring couple 
is only due to gravity, whereas in the second a torsional 
spring provides the entire restoring couple. 

The responses are evaluated approximately hy using 
the methods of harmonic balance [3] and multiple tjjne scale 
(MTS) [ ^2 ] . While using MTS method only the first order 
approximation is used and the results are c copared with the 
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first order approximate results obtained by harmonic balance 
method f It is shown that with the harmonic excitation, 
where the responses are also expected to be harmonic (and 
hence a priori knowledge of the solution is intuitively 
available), the res'ults obtained by the first order harmonic 
balance are accurate as compared to those obtained by using 
the same order approximation in MTS method. 

Solutions are also obtained by including higher 
oiSdfei' approximations in harmonic balance. Stability of 
the obtained solutions is investigated by perturbation 
technique and using Floquet’s theory [ 9] • Ihe results 
obtained by different order approximations are compared* 

It is shown that the first order approximation not only 
gives quantitatively inaccurate results in seme cases but 
,may even fail to predict the stability characteristics 
properly* The approximate results are also verified by 
numerical integrations of the equations of motion. 

For certain combinations of parameters the second 
order approximate results are found to be unstable. Only 
the existence of such combinations of parameters is sho^-m 
in this chapter. The study -of the responses under such 
situations is tmdertaken in the next chapter. 



2,2 SYSTM DESCRIPTION AND EQUATIONS OF MOTION 


Ilf 


Figure 2,1 shows the two degree - of - freedom 
system in which the primary consists of the mass M, the 
linear spring of stiffness and the viscous dashpot of 

damping coefficient c^ • The secondary system comprises 
of a simple pendilLum of length 1 under gravity field and 
hinged to mass M • The mass of the pendulum hob is m 
and k^ refers to the torsional spring resisting the 
angular motion of the pendiiLum. The pendulum motion is 
damped by a viscous damper of coefficient C2 • The primary 
mass is excited by a harmonic force F cos ojt » 

The equations of motion of the system are 

(M + m) X + X + k^ X - ml ( 0 sin 9 +0 cose ) 

= F cos “t ( 2 , 1 ) 

and ml 0 + 026 +^2® ** ml x) sine = 0, (2.2) 

•vdiere x represents the displacement of M ; 9 , the rota- 

tion of the pendulum, the dot implies the derivative with 
respect to time t and F and to are, respectively, the 
amplitude and firequency of the excitation. Equations ( 2 . 1 ) 
and (2.2) can be written as 

(1 + -§-) X + 2 i X - § 1 (e sine + e^ cose ) 


COS (I^t 


(2,3) 
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axLd 6 + 2 h + Q + ( ^2 ■" 1 ^ ® ~ (2.^) 

respectively, where 


/N 

^ is the natural frequency of the primary, 


■1 

*2 


2 M a 


is the daioping ratio of the primary. 


1 


y-^ is the linear natural frequency of the 
pendulum under gravity only. 


-4- is the natural frequency of the pen- 

■ (2.5) 

dulum with spring control only. 


(l3 ^ 


is the linear natural frequency of 
the pendulum under comhined control of 
spring torque and gravity 
c^ 


and C, 


2 ml^ to, 


is the damping ratio of the 


pendulum. 

When the pendulm does not oscillate, the system 
responds like a locked mass with the natural frequency as 




1 V M + m 

Equations (2.3) and (2.^) are written in the 
nondimen si onal form as 


( 2 . 6 ) 
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(1 + R) x"+ 2 




+X'-p‘^R(e sine + 

2 — 

0* COS0 ) =Pcost 


( 2 . 7 ) 


_2^ " ^ ^2 ^ P p -ij. 

and p 0 + -- e f 4. - ^ . - e + C ■:, — 7" p - - p^ x ) sinO 

/1+R 1+R 1+R-^ 

= 0 , (2.8) 

respsctively j where tiie variables are nondimensionalised as 


X = . and t = oit 


(2.9) 


The prime denotes the derivative with respect to the non- 
dimensional time t and the following nondimensionalisa- 
tions are used for the parameters; 


R 

q. 

and q 


m 


g is the mass ratio, p = 


0) 

IT 


P = 


F 


1 » 


(0/ 


“1 ’ 


11 = 


0 ). 


^^2 = 




(2,10) 


= / 


2 2 
^2 


At the excitation frequency equal to the locked mass natural 
frequency, p is given by 

p _ _J (2,11) 

/iTr 

When the system oscillates like a locked mass, 
i.e.,0 = 0 for t > 0, its motion is governed by 


(1 + R) p^ x^’ + 2C^ p Xj: + x^ = P cos t 
where Xj^ is the locked mass response nondimensionalised 




( 2 . 12 ) 
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with respect to 1 , The steady state solu‘:.ion of (2.12) is 

cos C t - ) , (2.13) 

where 

\ = — - ^ ' (2.lifa) 

y [1 - (1 +E)]2 + (2 p)2 

is the nondimensional amplitude of the locked mass and 

1 2 c. p 

$ = tan P — ^ 1 (2.1^b) 

. 1 - p (1 + R) 

2.3 APPROXIMATE SOLUTIONS BI l-iETHOD OF HAR40NIC BAL^CE 
AND STABILITY OF THE SOLUTIONS 

Equations (2.7) and (2.8) form a set of coupled 
nonlinear differential equations x^rhere the pendulum is 
parametrically excited. The solution is attempted in the 
neighbourhood of the primary unstable region of (2.8), i.e. > 

q 2 

with —5 — j- , when the static equilibrium posi- 
p^ (1 + R) ^ 

tion of the pendulum becomes unstable. Hereafter the 
static equilibrium position of the pendulum will be referred 
to as the zero solution of e . Furthermore* the range of 
excitation frequency considered is around tie natural 
frequency of the locked mass system, i.e. p /j + R Rf 1 

-1 

(from (2.11)), and q is taken as ^ . 
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As ttie excitation is harmonic, the steac3y state 
solution of tho primary mass is also assumed to he harmonic 
with, the same frequency oj ♦ The solution being sought 
near the primary unstable region of (2.8), the pendulum 
motion is taken as harmonic with frequency Thus the 

solutions are assumed to be of the form 

X = A cos ( t + ^ ( 2 . 15 ) 

and e = B cos ( 2 “ ^2 ^ (2.16) 

For accurate analysis higher harmonics should be 

considered in (2.15) and (2.16), in particular when the 

amplitudes are large. Substitution of (2*15) and (2*16) in 

3 9 

( 2 . 7 ) and ( 2 . 8 ) shows that some terms of order B-^, B etc. 
are neglected as they are associated with hi^er harmonics. 
But since we are investigating the existence and the stabi- 
lity of the harmonic responses with a harmonic excitation, 
tho analysis is carried out with the form given by (2.15 ) 
and ( 2 . 16 ). 

First the boundaries of the primary unstable region 
arc determined with the one term approximation given by 
( 2 . 16 ). 

2.3.1 Approximate Boundaries of the Primary 
Unstable Eegion 

For determining the primary unstable region, 
sin e is replaced by Q in (2.8) rendering if to the 


following form; 
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^ ^2 P ^ 

/I + R 


e' + [ 


1 + R 


2 -» - 
- p 25 ^ Je 


= 0, (2.17) 


where x-^ is obtained from (2.13) and (2.l4). 


Substitution of (2.16) in (2.1?) yields the first 
order approximate equation for the bounding frequencies of 
the primary unstable zone as 


2 

q2 p2 2 (?p p q) 

TTr ^ TTr~ 


If 

p 



If — » 


( 2 . 18 ) 


where is given by (2.14). 

For q = ^ , ( 2 , 18 ) gives the bounds of p for 
the primary unstable zone as 


I^boTJind 
where Q 


1 + R 


[ ( Q + 1 ) + + 2Q ] , 


(2.19) 


= — ( 


1 + R 


'2 ^ ^ 


+ ( 


1 + R 




The condition for these values of P^jound 


real, is 


P > 


2 c C 

2 

/I + R 


( 2 . 20 ) 


2 , 3*2 Approximate Steady State Solutions 

The steady state solutions of the equations of 
motion, (2,7) and (2.8), are assumed to be harmonic. 

The assumed forms are given by (2.15) and (2.16). aibsti- 
tutiing (2.15)..^d (2.16) into (2.7) and (2.8) and equating 
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the coefficients of the sine and the cosine terms j the 
following four equations in terras of the four unlmowns 
A, Bj and sire obtained: 


(v^ + ^2^ + 2 { 1 - (1 + R) } 

+ (2 p) v^] = ) 

+ ^2 , 


(2 21 ) 
( 2 . 22 ) 


tan 2 §. 


£ 1 - p^ (1 + R) ] V 2 - (2 p) 


V. 


( 2 . 23 ) 


[1 - p^ (1 + R)] + (2 p) V 2 + 


V- 


^2 

and tan (2 =7“ 3 

.2 ^ „2 


(2.2lf) 


B 

^ " 


B - 


where v 


1 


n^d.R) P^d^B) ^ 

(B) - (B) 


2 (B) 


?2 1 ® 

’yj' ^ ....ll ■■ lie . !£S .. I . I »i u ^ 

^ ” p </rT"R [ (B) + (B)]’ 

[ 2 (B) - (B)} 

+ B (B) + 2 J 2 <^B) +Ji^ (B)}] , 

00 ^ 

= £ 1 ~ (1 + R) ] + ( 2 p ) , 


and (B) is the Bessel’s coefficient of order n - arising 
out of expansions of the terms like sin (B cos ( ^ + § 2 ^ ^ 


. 
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Equation (2.21 ) is a polynomial in B, In the 
first order approximationj terms of order E onlj'" are 
retained in (B)’s. Equation (2.21) then reduces to a 
biquadratic equation in B uhich can be solved in closed 
form. In the second order approximation, terms up to B^ 
are included in (B)»s and terms up to B^ are retained 
in the third order approximation. In such cases the poly- 
nomial equation (2.21 ) can be solved only numerically. 

After evaluating B from (2.21 ), other unlmowns A, 
and are obtained from (2.22), (2.23) and (2.2^+), 
respectively. 

2.3*3 Stability of the Approximate Steady State Solution 

The stability of the assumed solutions, (2.1 5) 
and (2.16), is investigated by providing small perturbations 
to them. Thus the perturbed solutions are 

X = A cos (t + $^)+x^ (2.25) 

and e = B cos ( ^ + $2^ *^^2 * (2.26) 

where x^ and X2 refer to the small perturbations. 
Substituting these equations into (2.7) and (2.8) the 
perturbed equations are obtained as 

(1 + R) p^ + 2 P + x^ - p^ R [ f^ ( t ) x” 


+ f2 ( t ) x^ + f^ ( t ) X2] =0 (2.27) 
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and 


2 » 2 c p q ( 


/1 + R 2 
- p^ f, ( t ) 


+ R 


+ f 4 ( t ) ] 


X 


II _ 


= 0 


(2.28) 


respectively, 
( t ) 

C t ) 

. f.^ ( t ) 

and f ^_ ( t ) 


xjhere 

= sin [B cos ( ^ + $2^ 3 > 

X 

= - B sin ( 2 ^2 ^ ^ ® 

B t 

= - jj- cos ( 2 ■*■ ^2 ^ t ® 
2 

- ^ sin^ ( ^ ^ t ^ 



( I +«2)] 
( 2 $2^ ^ 
( I +* 2)1 


COS [ B cos ( 


* 2 )] 


( 2 . 


Equations (2,27) and ( 2 . 28 ) are linear equations 
with periodic coefficients. The functions f^^ ( t )*s in 
( 2 . 29 ) are computed with the values of A, B, and §2 
obtained by the approximate results. The order of the 
terms of B, in expansions of f^ ( t )’s, are kept same as 
those in (2,21 ), Floquet's theoiy[ 9 ] is used to investi- 
gate the stability of (2.27) and (2.28), Qiiese equations 
are expressed as four first-order differential equations 
in the following form: 


29) 



•where {w) 


T 

= {x^ x^ x^ x^ > 

f ( t ) is a periodic matrix with period 4 is « There- 
fore, if[WCt)j is a fundamental matrix solution, 
[W(t+ifis )J is also a fundamental matrix solution. 

Let these he related hy 

[ W ( t + )] = [ S] [W ( t ) ] C2.?.0) 

To generate [ S] numerically, (2.30) is integrated with 
initial conditions [W (O)] = [l] for one period of 
oscillation, where [ I ] is the identity matrix* Then 

[ S] = [ W (IfTt )] (2.31) 

The eigenvalues, , of [ S ] are the roots of 

I [ S] - X[ I]( = 0 (2*32) 

The steady state solutions, (2.15) and (2.16), are stable 
if |X^| < 1 for i = 1 to 

Though numerically this process is involved, 
this is justified since in seme cases the higher order 
approximations are required to predict the true stability 


behaviour. 
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2«3»^ First Order Approximation by Harmonic Balance 

In the first order approximation, terms of order 
B only are retained in (B)‘s in (2.21 ) to C2.2if). 
These equations then bee one 


^ 2' 


(p"^ R + 16 R B^ + 6^ [ "'^6 ** ^ 3 = 

2 (2.33) 

P It- 

A = + i ’ ^hI r, ]) 


,2 4 


(2.34) 


[1-/ (1.H)]^-2SpC^-T^) 


tan 2 §2 


P ^ / 2 2 

and tan (2 S ^ ^ -j ^ " ' ^" ■ " y ^ ^ - Y + 1 ^ » 


respectively, where 

2 2 

■Vi- = [ 1 (1 + R) ] ( ^ - TTr ^ 


(2 p) 


^2 P ^ 

/TTr 


(2.35) 


(2*36) 


(2.37) 


. p2 ci2 ,2 ( ?2 P 'i)' 

and vg = ( -q: - tpjTR ) + 1 +R 


Equation (2.33) gives two roots of B as 


[ - ± - ^4. ^6 3 (2.38) 


2 „if 


p^ R 
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Equation (2,l8)) giving the bounding frequencies of the 
primary unstable region, can also be obtained by putting 
B = 0 in (2.33)* It can be shoim that the inequality 

^ holds within the primary unstable 
region when (2,38) yields only one real root of B as 

Upto certain values of p , on either side of the primary 
imstable region, (2.38) yields two real roots of B if 

2 ^ 

■^5 + ^6 > ^ (2.4-0) 


Moreover, in these ranges of frequency the zero solution of 
0 is stable. Such frequency zones are termed as 'overhangs*. 
The inner boundaries of these overhangs are simply the boun- 
daries of primary unstable region. The outer bounding 
frequencies for the overhang regions are obtained from (2<4-0) 
as 


1 + R 


where 


= 1 + 


[ Q.J + 




1 ] , 


( 2 . 41 ) 


— J + c ^ ) 

1 + R ^ 2 ^ 


2 ( 


It is obvious from (2.4-1) that if = ^2 ~ 

the overhang regions cover the entire frequency range (of 
course excluding the primaiy unstable region). The first 



order approximation yields same overhang regions for the 
gravity controlled and the spring controlled pendulums j 
which obviously is not true* It is shovm later that the 
first order approximation gives only a rough estimate of 
the overhang regiens. 


2 A MULTIPLE TII4E SCALE METHOD : FIRST ORDER APPROKIMATIOK 

To deteimine the s’oitability of harmonic balance 
method for the present problem, results arc also obtained 
by anotlicr pov^erful approximate method namely multiple time 
scale method [42] . The results obtained by these two 
methods are compared carO-y for first order approximations* 
To tliis end (2.?) and (2.8) are written as 
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X ( t ) = e ( Tq , T^) + e ^ (Tq > (2A5) 

and e ( t ) = e ( Tq , T^) + (Tj, , T^) , (2.lf6) 

•idle re 

Tq = t j = et (2A7) 


To order the forcing and damping terms so that "they appear 
in the same per'turhation equation as the nonlinear terms, 
\TQ introduce 


and 


’2 

P 


= e P- 2 

2 — 

= 2 £ P 


% 


> (2.if8) 


Further, the solutions are sou^t for to near to^ and for 

i 

the frequency ratio q = ^ > when 


p = 1 + ecj 

and q = ^ 


(2.if9) 


Substituting (2.4-5) - (2.4-8) in (2.4-2) and (2.4-3) 
and equating coefficients of like powers of € 

X. + X. = 0 ( 2 . 50 ) 

o 1 1 

Dq + ^ = 0 (2.51) 

Dq + ^2 = - 2 cy - 2 - 2 x^ 

+ 2 P cos t 


(2.52) 
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ana = _ 2a y^ . 2 y^ - y^ 

+ (Dq , (2,^3) 

^hore D^ = -|^ , d 2 = and D, 5 ^ . 

O ^ *c , 1 

The solutions of (2.50) and (2.51 ) are ey.pressed as 

= A (T^ ) exp (i T^) + A* (T^ ) exp (- i T^) 

i T i T (2.^) 

and y^ = B (T^) exp ( ) + B* (T^) exp ( - — 

respectively, ^^rhere A and B are the canplex conjugates 
of A and B. Substitution of (2.5^) in (2*52) and (2.53) 
yields the solvability condition as 

•D 2 -jS _ 

2crA-2iD^A-2i l^^A - + P = 0 . 

rt -D B _ (2.55) 

and ^ - i B - i — ^ - A B = 0 

Introducing the polar coordinates 

A (T J = I a (T^) exp [i ) : 

(2.56) 

and B (T^ ) = -^ b (T^ ) exp [i Y ^ )^ 

in ( 2 . 55 ) and sepj^rating the real and the imaginary parts 

a = - a + % sin ( - 2 Yg) _ p sin ' 

"R — 

aD^Y^=*- (S d- + cos ( y j - 2 Y 2 ) *" B cos Y ^ 

b = - - 1^ sin ( Y^ - 2 T 2 ^ ' 

and ^2“ - *^+2 i t ^ ^ 2 t 
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For steady state solutions 


a = D, = D, b = I>1 1^2 = 0 


( 2 . 58 ) 


and hence (2,57) yield 

2 
^2 


o2 2 2 

O' = CJ + (J. ^ 


( 2 . 59 ) 


and i i a'^ \x ^Vl 2} ± 


r. 


p - <y ( + iig) J 


( 2 . 60 ) 

To check the stability of the steady state solu- 
tion the coordinates are perturbed as 

a = 3 .^+ 3 .^ , 


b = b^ + b^ , 


= Y^o ^ 


and ^ 2 ~ 20 ^ 2”! ^ 


> (2.61) 


where suffix 0 implies the steady state "values. Subs- 
tituting (2,61) into (2.57), expanding for small a.j, b^ 
etc. and keeping only the linear terms in b^ etc. the 

following equations are obtained: 

■)3 

a^ = - ^ ^ ^10 " ^^ 20^^^1 

b2 

+ ■# ^ "^10 - 2 Y ^ q ) 


- P cos Y^q ] Y-ii ■ 

b2 

-Cr^ -i <=“= ^ ^0 - 2 


(2.62a) 
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r !R. ^ ^ 

Y-i = - L cos ( - 2 Y20) - P cos Y-jo 

b 

+ iri;=°« < 1^10 -? 1^20 

r R ^o ~ - 

1 + R 8 a^ ^ ^10 ” ^ ^20^ " ^ ^10 ;l^11 

^2 

twR iri, =“ ( 1^10 - 2 1^20^^21 

r -, a b 

Dl = -[^ sin ( Y^o - 2 -[-^^ cos ( ^ 

+ [^0 ^o ^ ■'^10 " ^ '*^20^^ ^21 (2.62c) 

^1 ^2 = t 2 ^ ^10 “ ^ “^20^^ ^1 ‘■^^“2 ^ Y-^0 " ^ ^20^^ ’'^11 

+[a^ sin ( Y-jo “ ^"'''20^^21 (2.62d) 


The stability of the steady staiie motion depends 
on the eigenvalues of the coefficient matrix on the right - 
hand side of (2.62). These eigenvalues, are the roots 

of the polynomial 

+ d^ X^ + d^ x^ + d^x + dj^ = 0 , (2.63) 

where d^ = 2 + ^-2 j 

' ^2 

d2 = -1 ( + 2 + a + ^ ^ ~ 2“ ^ > 

o a^ 

b2 

43 = H I + (5 1*2 + 3 + 2) ] 

0 a^ 

o 2 

+ (|>^ + cr ) P2 



As , d^ and d^ are positive, the necessaiy condition of 
stability is 

dj^ > 0 (2.6V) 

Tliis implies that in (2.60) the plus sign before the quan- 
tity under the square root gives the stable root and the 
minus sign gives the unstable root. 

2.5 KESULTS Al'ID DISCUSSIONS 

The results are obtained with R = 0.2, = 0.02, 

= 0 . 05 > q = 0.5 and varying the other parameters q. jq 25 ? 
and p in a suitable manner. Initially the first order 
approximate results, obtained by HTS and harmonic balance, 
are discussed. Then the higher order approximate results, 
obtained by harmonic balance only, are presented. The 
approximate solutions are comparod with the niunerically 
integrated results. For this, (2.7) and (2.8) are numeri- 
cally integrated using Gill's modification of tiie classical 
Runge - Kutta method [ Wj . Witli the period of forcing 

71 . ^ . 

function as 2^ , a step size of is found to give 

good results. All the results in this section are obtained 
with this step size. Inside the primaiy unstable region, 
the initial conditions are taken as x (O) = x* (O) = 

9'(0) = 0 and e(0) = 0 . 05 . In the overhang regions 
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the integrations are performed Toy clianging the frequency, 
p, in steps when the solution has reached steady state. 

2,5*1 First Order Approximate Results 

In the first order approximation, the amplitudes 
of X and Q are calculated by harmonic balance using 
(2.3^) and (2.38). The same are also obtained by using 
MTS \hLth ( 2 , 59 ) and (2.60). It is obvious that so far as 
the first order approximation is concerned the results 
depend only on the value of q. and not on the Individual 
values of and fact the first order approximation 

(by both the methods) essentially gives the solution of the 
approximate equations of motion. These approximate equations, 
obtained by putting sin e = e and cose =1 in ( 2 .?) 
and ( 2 , 8 ), are 

(1 + R) p^ x* + 2 p X* + X - p^ R (0 0 + 0 ^) 

= P cos t ( 2 . 65 ) 

and 9 ’ ^ 0 * + ( x" ) 0 =0 ( 2 * 66 ) 

Equations of this form with quadratic nonlinearities have 
been studied in detail, for exajxsple, in references [9j 26 ] . 

Figures 2.2 and 2.3 shov/ the amplitudes A and B, 
respectively. The results obtained by MTS are shown by 
the curves 1 - 1 and thoso by harmonic balance are given by 








FIG* f.3 VABIATION OP AMPLITDBE CF PENDULUl-l WITH 
EXCITATION FBBQDiKCI ; FIRST ORDER 
approximate results, P =0,02 

4 
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the curves 2-2, The fin: lines show the amplitudes x/ith 

frequencies lying inside the primary unstable region and the 

dotted linos indicate the stable solutions ia the overhang 

regions. In MTS the stability condition is given by 

(2,6^), which states that the upper branch in the overhang 

(shown by dotted lines in Fig, 2*3) is stable and the loxjer 

branch (not shown in this figure) is unstable. It is obvious 

a- 

from (2,59) that no sepj^rate lower branch exists in Fig, 2,2* 
The stability of the solutions obtained by harmonic balance 
is chocked with (2.32) and the resiiLts are same as those 
obtained by MTS. The responses obtained by integration of 
(2.65) and (2.66) are also found to be harmonic and the 
amplitudes are indicated by the circled points in Figs. 2,2 
and 2*3» 

Both MTS and harmonic balance yield same results 
at p = 0,913 ( p = 1.0). While the response curve 'obtained 
frem MTS is symmetrical about p = 0.913} that by harmonic 
balance is not so. The results by these two methods start 
deviating as p is varied ax'jay from 0*913 on either side. 
The harmonic balance results are seen to be closer to the 
integrated results. This is due to the fact that in 
hajTOonic balance no approximation is made on the frequency 
range, whereas MTS gives resxiLts in a narrovx frequency 
band as given by (2.^9)* 3h fact, the same approximation 
on the range of frequency makes (2,^) and (2.38) same as 
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(2.59) and (2,60), respectively. 

If the excitation parameter is increased the 
amplitudes will also he large and the first order approxima- 
tion will no more be valid. Moreover, it is obvious tiiat 
with large amplitudes the nature of the pcnd-uLum motion (and 
hence the motion of the primary mass) will depend on the type 
of the restoring couple it is subjected to - a feature which 
can not be revealed by the first order approximation. Since 
the frequency range of interest are quite x^ride, the subse- 
quent higher order approximations are carried out by the 
method of harmonic balance only. 

2,5.2 Higher Order Approximations by Harmonic Balance 

In the first order approximation only one value 
of A is obtained in the overhang regions as given by (2*3^) 
and (2,59). Further, the first order approximation renders 
A independent of uho excitation parameter, P, as also noted 
in references [9, 20], In the former reference this was 
referred to as saturation phenomenon. This saturation pheno- 
menon is the consequence of the first order approximation and 
not a system characteristic, Eiis is revealed from Fig, 2,4- 
which shows the first order approximate results obtained by 
( 2 . 34 ) and ( 2 . 59 ) and the integrated results of (2,65) and 
(2.66) at the locked mass natural frequency (p = O.913). 

All these results arc same (shovm by tiie dotted line in 
Fig. 2.4) and show the saturation phonemonon for P > 0,0018. 
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The chain dotted line shows the amplitude of the locked mass 

2 ?r 5" 

response. For P < 0.0018 (l.o.,P < - vide (2.20)) 

, , . , . ... / 1 + E 

the static equilibrium position of the pendulum is stable. 

Only for P > 0.0018 the pendulum motion affects the primaiy 

response. Also included in this figure are the results 

obtained by the second and 1iie third order approximations 

(using (2.21) and (2.22)) and by integration of the original 

equations of motion (i.e. (2.7) and (2.8)), It can be seen 

that in these cases the amplitude. A, is no longer independent 

of P and varies almost linearly with P for small values 

of the excitation. 

Figures 2.5 and 2.6 show the variations of the ; 

amplitudes, B and A, respectively, with the nondimensional forcin 
frequency, p. The curve 1-1 indicates the steady state 
amplitudes obtained by numerical integration of (2.7) and 
(2.8) with q.j “ 2 * ^2 ~ ^ (i.e., the gravity controlled 

pendulum) and the cuive 2-2 shows the integrated results 

1 ! 
with q.^ =0, ^2 ~ 2 spring controlled 

pendulum). The circled and iiie squared points indicate the | 

stable amplitudes obtained by the first order and the second ^ 

order approximations, respectively. The primary unstable | 

region obtained by (2.19) is represented by PTIR, where the 

approximate solutions yield only one stable steady state. , 

The dotted lines indicate the stable amplitudes in the 

I 

ovcrhsing rBgionsj th© zGro solution of 0 is also .static • 
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2-2 ; -i -! 

r- INTEGRATED RESULTS i4J,.L:i 

O RRST ORDER APPROXIMATION , . ;. ] J ! 

,a SECONDORDER APPROXIMATION . T. 

LOCKED MASS RESPONSE : i 





The 3 imps at the ends of the overhangs are showi hy arrows* 
The unstable solutions in the overhangs are not shown. The 
extent of the overhang region obviously depends on the type 
of the pendulum's controlling torque. It is evident that 
both the amplitudes (A and B) are obtained very accurately 
by using the second order approximation. The integrated 
responses are also fomd to be completely harmonic. 

If the excitation parameter, P, is increased 
beyond a certain value, the steady state solutions are seen 
to deviate from being periodic x^ith a predoninant first 
harmonic as assumed in (2,l5) and (2.l6). This critical 
value of P depends on the values of other parameters. In 
such cases the characteristics of the solutions are different 
depending on whether the restoring couple csn the pendulum is 
due to gravity (i,e,,q2 = O) or due to the torsional spring 
(i,e.,q^ = O), These two cases are separately discussed 
^>€low, 

(1) Case I : q^ “ 2 * I2 ~ ^ 

The approximate results with P = O.O^ are presen- 
ted in Figs. 2.7 and 2.8, where the numerically integrated 
results are also shown by firm lines. The first, second and 
third order approximate stable solutions are shown by the 
circled, squared and triangled points, respectively. The 
second and the third order approximate results have negligi- 
ble difference up to p = 0.95 a-O-d match well with the 










integrated results* Numerical integrations also reveal that 
with p near 0,955 ’v^en 8 is given a slight distiirbance 
(say in the form of an initial condition 6(0) = 0,05) the 
pendulimi makes seme revolutions during the transient state 
and ultimately settles down to a periodic steady st&te. 

These number of revolutions are seen to depend strongly on 
the initial conditions. 

As p is increased beyond 0*95 the difference 
between the results obtained by the second and the third 
order approximation becomes appreciable, signifying the 
absence of the dominant first harmonic in the steadj'' state 
responses. The stability analysis using (2.27) and (2.28) 
reveals that these solutions are unstable for p^ “1 *0 
thus implying that the assumed form of solution may not 
hold good even for such lo\7 values of the excitation para- 
meter, P. In such cases the numerically integrated responses 
are also found to be nonperiodic. These responses show a 
complex wave form and are extremely sensitive to the initial 
conditions. A continuous and irregular transfer of energy 
between the two modes is apparent from these complex wave 
forms ^ 

It can be noted here that for p >_ 1.0 (upto 1.30) 

the first order approximation by MTS and harmonic balance 
yield stable solutions with nonzero value of B. In such 
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cases the integration of (2.7) and (2,8) yields nonperiodic 
responses for 1,0 £ p< 1.1 (after which the zero solution 
of 9 is strongly stable), whereas, that of (2.65) and (2,66) 
yields periodic responses (with nonzero araplitude of 0 ) 
for 1,0 ^p< 1.30. Further, the first order approximation 

by harmonic balance shows that the approximate soluticais are 
unstable for p< 0,78 which is again verified by integra- 
tion of (2.65) and (2.66), Hence it may be observed that the 
first order approximation reveals the behaviour of the approxi- 
mate equations, (2,65) and (2.66), rather than that of the 
original equations of motion. 

The responses of the system with excitation para- 
meter P ^ 0,05, where the responses are nonperiodic over 
certain frequency range are studied in detail in the follow- 
ing chapters. 

(2) Case II: q^ = 0, 0-2 = 2 

Figure 2,9 shows the amplitude B for this case 
with the values of the other parameters kept same as in 
Fig, 2*7 » The approximate results are seen to be quite 
accurate for the frequencies well inside the primary unstable 
region (PUR), In the present case, numerically integrated 
results show periodic steady state responses at all frequen- 
cies, It should be noted that, unlike in the previous case, 
complete revolution of the pendulum is prevented by the 













spring* This in fact renders the steacJy state solutions to 
he periodic in the entire frequency range* 

ihe approximate solutions differ from the integra- 
ted results near the two peaks. Integration reveals the 
existence of two steady state solutions at the left peak 
(p =: 0*83) one where the resxjonse is almost harmonic 
(shown hy point A in Fig, 2*9) and the other where tlie res- 
ponse is amplitude modulated (shown hy point B), In the latr 
tor ’ case (point B), x response is found to have small 
modulation with period 4 75 and 0 - response is found to 
include some even harmonics as v^ell. At p = O.83, the 
response time histories for x and e are plotted in Fig. 
2*10. Case (a) shows the steady state responses obtained 
hy decreasing p in steps from p = 0.84 and case (h) shows 
those obtained hy increasing p from p = 0,82. The pen- 
dulum response has higher harmonics even when p < 0,83* 

Near the right peak also integration reveals the 
amplitude modulated response of x and e - response has 
other terms present beside the first harmonic. In such 
cases the amplitude hy integration (Fig, 2.9) represents 
the maximum value during one period of oscillation. Because 
of higher harmonics, tlie form of solution, assumed in (2.15) 
and (2#l6), is not correct and tiiis explains the difference 
between the approximate and the integrated results near the 
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two peaks* At seme frequencies near these peaks (e.g* 
p = 0«81 and 1*10) no stable solution by liaimonic balance 
is obtained. Only the amplitudes of the stable solutions 
by the approximate method are shovm in Fig, 2*9* 

Lastly, it may be noted that at p = 0.83 (the 
critical frequency) the bifurcation of the stable harmonic 
motion to one unstable harmonic motion and another stable 
amplitude modulated motion is associated \;ith the jumps in 
the amplitude values, as shown in Figs. 2,9 and 2,10. 



CHAPTER 3 


RESPONSE WITH LARGE EXCITATION 


3.1 INTRODUCTION 


In this chapter the holiaviour of the system, 
already considered in Chapter 2, is investigated with large 
values of the excitation parameter P. The typical range 
of P studied in this chapter is 0,05 < P X 0*15* Here 
only the case of the gravity controlled pendulum is consi- 
dered, i.e., = 0 and q = q^ in (2,8). The equa- 

tions of motion, thus, are same as (2.7) and (2.8) and are 
reproduced below; 


(1 + R) 


2 

p X 


.1 — 


+ 2 +X-P R(e sine -f 


9^ cose ) 


p 


2 


9 + 


2 q p 

/TTr 


e + ( 


= P cos t 

Op- 2 — ” s 

- p X ) sin e 


(3.1a) 


= 0 


13 . 113 ) 


The linear naixiraT ■cai^quency of tne pendulum is 

t/ccE-en as half of that of the locked mass system, so that 

1 ' 1 
q = All the resTolts in this chapter are for q = p . 

First the boundaries of the unstable regions of 
the zero solution of (3. It)) are discussed. Inside the 
primary unstable region, steady state solutions are investi- 
gated using the method of harmonic balance. The stability 
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characteristics of these assumed harmonic solutions are then 
checked hy Floquet's theory. Inside the primary unstable 
zone two different types of regions are discussed - one 
where the assumed haimonic solutions are stable and the other, 
where no stable harmonic solution exists. The boundaries 
of these regions are also obtained. 

These results are then verified by numerical inte- 
gration of (3*1). First, the steady state amplitudes found 
by approximate method and numerical integration, are compared 
for the cases where (3*1) has a steady state with x and Q 
oscillating harmonically. For the regions where the approxi- 
mate method does not yield any stable harmonic solution, 
only the numerically integrated results are presented. Such 
regions are classified in two categories - (i) where there 
may still be one or more periodic steady states depending 
on initial conditions, and Cii) where no such periodic state 
is obtained. For the later case, some physical justifica- 
tions are given to analyse the results further, even though 
the integrated results may not converge. Lastly some inte- 
grated results are presented briefly for sufficiently large 
excitation. 
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3\2 BOUNDARIES OF THE UNSTABLE REGIONS OF 
THE ZERO SOLUTION CF 0 


For studying the boundaries of the unstable regions 
of the zero solution of e , sin 9 is replaced, by 9 in 
(3.1 b) rendering it to the following form: 


P 


2 


It 

e 


2 p q 

+ — 9 

/ 1 + R 


1 + R 


2 — ” \ 
P^ )6 


0 ( 3 . 2 ) 


Here Xr is the locked mass response and is given by (2.13) 
Ju 

and (2.1^). 


The boundaries of the primary instability of (3*2)} 
with the first order approximation, can be found from (2.19). 
However, for large values of P, the one tern solution assumed 
in (2»l6) does not yield correct bounds and more number of 
terms in the solution are to be considered, i.e., with q - 2 > 
0 should be assumed as 


n — S 




[ 




it 


Instead of going through this process, (3»2) is 
transformed to the standard form of Mathieu equation and the 
boundaries of instabilities are determined with the help of 
the standard tables. Substitution of (2.13) in (3.2) yields 


2 P q 


P 0 


0 


+ R 


+ [ \ ^ t = 0 


o 

C3.4) 



5 ^ 

To convert (3.4) into the standard form of Maiiiieu equation, 
the following transformations are used: 


t - = (2 u - ft) (3.5a) 

2 q u 

and e = exp ( ) (3.5b) 

p /I t R 

Substitution of (3*5) in (3.4) results in 


^ ( c . 2P 

d u 


cos 2 u) ^ G j 


where 


a = 


and 

P = 


4 q^ (1 - ) 

p^ (1 + R) 


^ v/[l - p^ (1 + R)] ^ + (2 p)^ 


(3.6) 


(3.7) 


Equation (3.6) is in the standard form of Mathieu 
equation [l j . The unstable solutions of (3*6) are bounded 
by the periodic solutions of periods te and 2^ alternately. 
The primary unstable region is bounded by two periodic solu- 
tions of period 2 — one of them being odd and the other 

being even. The odd solution is represented by b^ (se^ ) 
and is given by a sine series. The even solution is repre- 
sented by a^ (ce^ ) and is given by a cosine series [ 1 ] . 
Similarly, the secondary unstable region is bounded by two 
solutions of period ^ , one of them is odd, represented by 
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other being even, represented by ( 062 )* 

The boundaries of instabilities of ( 3 . 6 ) are 
determined directly from the tables [45] • These tables are 
given for the following cqiation; 


d u" 


( cos'^ u ) = 0 


(3.8) 


Therefore, noting the values of and from those tables, 
the boundaries of (3 .6) are determined with 


and 


P = 


"TT 



(3.9) 


From ( 3*5 b) it is obvious that because of the damp- 
ing in the pendulum, ?2 > "^^e regions of instabilities of 9 
will be slightly narrower than those of , It was shown 
in ( 2 . 20 ) that to cause instability with <1 = -5 , ? should 

2 C. Cp 

be greater than — . Since the values of 5 ^ and 

/1 + R 12 

are quite small, the boundaries for 8 a.nd i{) will be almost 
same, except for very small P. 


The bcamdaries of the primaiy instability region 
of (3.6) are shown in Fig, 3.1. The unstable regicn is 
bounded by b^ (se^) and (ce^). For a given excitation 
parameter, P, and forcing frequency, p, the correspanding 
values of a and P are calcialated from (3.7). 




F3B. 3*1 REGION OP PRIMARY INSTABILITY OP MATHIBIJ BQUATICW 











For a certain value of P, the range of p over 
which the zero solution of 9 is unstable, is found hy 
plotting the locked mass response curve in the a - 0 plane 
in Fig. 3*1 • The points of intersection of this curve and 
the instab ilitj^ boundary give the bounds on a . Ihe bounds 
on p are then determined using the relation 


2 ^ (1 - 
P = a ( 1 + R) 


(3.10) 


For easier inspection, it is desirable to plot the 
boundaries of instabilities in P — p planev To this end, 

P is expressed in terms of p as 


p = -f Jll - p^ (1 + R) ] ^ + (2 C. p)^ (3*11) 

For each point ( oc , P ) on b-j (so^ ) etc., the correspondino 
point in P - p plane is determined uslig (3.10) and (3*11 )• 
These bomdaries are plotted in P - p plane sno^vn in Fig. 3*2, 
where b^ (so^) etc. denote the correspondjng boundaries in 
a - p plane. From (3.5a) it is clear that b^ (se^ ) and 
(cG.) have period hn , bp ( 362 ) and &2 having 

a ‘period 2% • 


3,3 APPROXH-IATE SOLUTiaiS IN THE 

UNSTABLE REGION AND THEIR STAB IL HI 


The method of harmonic balance, explained in 
section 2.3, is used to obtali approxlnate solutions of 
(3.1) in the primary unstable region. It was observed in 
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SGction 2*5*2 that for P = 0*05? the first order approxi- 

mation does not yield correct results, quantitatively as well 
as ' qualitatively • Therefore hi^er order approximations ai-e 
considered in this chapter by including higher power tonus 
in the expressions of Bessel's coefficients, (B)'s, used 
in (2.21 ) to (2.2^)* Since the amplitude of 8 is expected 
to he quite 'ligh for largo P, the tliird and the fourth order 
approximations are used, i.e., terms upto and B'^ , 
respectively, are considered in (2*21 ) to (2.2^). The solu- 
tion, sought for large values of P, is still harmonic and is 
of the form given by (2.15) and (2.16), 

The stability of these solutions is investigated 
by PloquGt's theory after obtaining the perturbed equations 
(2,30). 

3, if SOLUTIONS THROUGH HUJIERIGAL INTEGRA.! ICNS 

Equations of motion (3*1 ) are integrated numerically 
using Gill’s modification of Rungc-Kutta method [ if6, Mf ] . 

For the cases where the solutions roach periodic steady state, 
a step size of integration equal to is found satisfac- 

tory, In some cases the step size is further reduced to 
^ to check the validity of the results. "When the results 
do not converge with decreasing step sizes, integrations are 
also performed using Hamming’s modified predictor - corrector 
method [46, 4-7] . For bettor accuracies, the double precision 
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ccmputations are also carried out. The upper error hound 

-3 

in the predictor - corrector scheme is varied from 10 to 
10 The initial step size is tahen as ^ • If more 

than ton bisections of the initial increment are necessary 
to get satisfactory accur*acy, the integration procedure is 
terminated. 

The results presented in this chapter are obtained 
by Gill’s method. Special references are made to the predic- 
tor - corrector method wherever used. 

3.5 RESULTS MD DISCUSS lOHS 

The characteristics of the solutions (discussed in 
this section) depend on the combinations of amplitude, P, 
and frequency, p, of excitation. These combinations in turn 
depend on the values of other parameters* All the results 
presejnted in this section are for the following values of 
the parameters: R = 0,2, 0.02, ^2 ” 

q = 0,5* 

3 . 5.1 Boundaries of Instabilities of the Zero Solution of 0 

The primary instability boundaries given by 

a icej and b^ (se, ) are plotted in Eig. 3*1 in a - p 
11 11 

plane. The locked mass response curves, hereafter referred 
to as Lt-IRC, are also plotted by dotted lines in the same 
figure for v8.rious values of P. For lower values of P, 
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(say P - 0,05), LMBC intersects th.e primary boundaries 

only at two points. These points of intersections determine 

the rounds on 0 . as oc = 0.73 and oc = 'I •3^5- 

Fresn these values, the hounds on p are calculated using 
(3.10) as = 0.78 and = 1.067. 

For a large value of P (e.g. P = 0,2), however, 
it can he seen from Fig, 3,1 that LURC intersects a.^ (ce.^) 
at three points. Near a= 1? L^'IRC passes on to the sccondfiry 
region of instability. This is more clearly seen in 
Fig. 3*2 where the boundaries are plotted in P - p plane. 

In tliis figure b^ (so.|) and (cc.j) bound the primary 
unstable region, and and a2 ^^^2^ 

boundaries of secondary unstable region. For P > 0,12, 

near p = 0,913 ( p = 1*0) the secondary unstable region 
starts appearing and this region expands till P ^ 0,36, 

For P> 0.37, the third unstable region appears which is 
not shown in the figure. 

In this figure the boundary of primary Instability 
obtained by the first approximation, i,e., from (2,19), is 
also shOT/jn by the circled points. It is clear that the one 
term approximation of (2,16) does not give good results for 
P > 0.10. 



3*5*2 Results by Approximate Method 
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Ihe approximate results) presented in this section, 
are obtained by using the fourth order approximation. These 
results do not differ much from those obtained by using the 
third order approximation, 

^ig* 3*3 shows the enlarged view of the boundaries 
of instabilities, already given in Fig, 3,2, for P upto 0.15. 
Because of dampings and primary boundary starts 

from a nonzero value of P (0.0018), as dictated by (2.20), 
The boundary of the primary instability is shown by the firm 
line. For sufficiently large P and certain range of p, 
the approximate solutions are found to be unstable. The 
zones where the approximate solutions are unstable are found 
to be seperated by the zone of stable solutions. These 
zones of stable and unstable solutions are shown in Fig. 3*3* 
The t&>stable zones (regions II, III and IV) are shown by 
hatched areas. The stable zone is indicated as region I, 

In region IV the nonzero solutions are found to be unstable 
whereas the zero solution is stable. The areas between the 
firm line and the chain dotted lines show the overhangs 
where the zero solution is stable. The nonzero solutions 
in the overhangs may or may not ho otable depending on 
whother the representative point lies in region I or not. 

As the entire region IV is quite far away frcsn the primary 
boundary, the static equilibrium position of the pendulum 
is strongly stable and the system is expected to behave 











like a locked mass systsn. li the discussions to follow, 
the nature of the solution in region P/ is talccn as that 
of the locked mass and is not considered any further. 

It is obvious from Pig* 3*3 that as the frequency 
p, is varied, the representative point passes tlirough 
various regions depending on the value of P. As an example 
we consider the behaviour of the solutions for P = O.07, 

For 0.685 j< p < 0.84- the soXutions are stable. As the 

frequency is increased beyond C.84-, region II is entered 
and the solutions are unstable. After that for a narrow 
band of frequencies, 0,89 < p O.96, the representative 
point is again in region I and the solutions are stable. 

For further increase in p, region III is encountered 
rendering the' solutions unstable. Another interesting 
feature to be noted from this figure is that the right hand 
side overhang diminishes with increasing P and ultimately 
disappears for P > 0.08, The left hand side overhang also 
becomes narrower for largo P and in tliis overhang the 
solutions are stable only in a narrov; range of p. 

Finally, with P > 0.12, n^ar p = 0.913 Cp = 1-0) 
secondary region of instability starts appearing. Since 
the secondary region of instability is bounded by the solu- 
tions of period 2 , the assumed form of solutions, given 

by (2.15) and (2.16), is not correct. However, this regicn 
of secondary instability is not analysed by the appraximate 


method. 
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The discussions presented so far dealt v/ith the 
qualitative aspects of the solutions with special reference 
to their nature of stability • We now turn our attejntion 
towards some quantitative aspects of the responses# To this 
endj Figs# 3#^ and 3*5 show the amplitudes of the primary 
and the secondary, respectively, for different values of 
the excitation parameter. The amplitudes, corre spending to 
the stable solutions, are joined by firm lines, whereas, 
those corresponding to the unstable solutions are encircled 
or squared and are joined by dotted lines. The squared and 
the circled points refer to the solutions lying In regions 
II and III of Fig, 3#3} respectively. The overhangs are 
included in these figures but not indicated specifically# 
Termination of a response curve with firm line indicates 
the end of the overhang, after vjhich the zero solution of 
0 is strongly stable. On the other hand, if a response 
curve ends with a dotted line then it moans that the entire 
range of p, for which the zero solution of 9 is unstable, 
is not covered in the figures. 

It is seen from these figures that for P > 0.05, 
the amplitudes sharply Increase, at certain frequencies# 

In fact, these are the frequencies at which region III 
(Fig, 3,3) is approaheed. Moreover, because of tire ^arp 
increase in the amplitudes of 6 there is an appreciable 
difference in the successive approximations while solving 
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(2*21) It can te remembered that for P = 0.0^, 
this point was discussed in section 2.5.2. Hear the frequen- 
cies at which region III is approached, there may be more 
than two solutions of (2.21). As an example, the niimber of 
solutions of ( 2 * 21 ) and their stability characteristics 
with P = 0,07 at frequencies near region III are as follows: 

(i) for 0,89 <_p 0.95 there is only one stable solution, 

(ii) at p =: 0,96 there are three soluticns, only one of them 

being stable, (iii) at p = 0*97 there are three unstable 
solutions, and (iw) for p > 0,98 only one unstable solu- 

tion is obtained. The values of the roots of (2.21) and 
corresponding values of A for above mentioned cases are 
given in Table 3*1 • For the cases v/here only unstable 
solutions are obtained and the number of such soluticns 
exceeds one, that solution which gives a smooth curve for 
the amplitudes is shown in Figs. 3*^ and 3*5* Fhe ampli- 
tudes of these unstable solutions increase sharply 

and then start decreasing as the frequency p, is further 
Increased. This decreasing trend, however, is not shown 
in these figures. The complete amplitude curves are shown 
only for P = 0,03. At this value of excitation the 

steady state approximate solutions are stable for all values 
of p. From the amplitude curves for P = 0.03 it can be 
seen that both A and B increase slowly after p = 0.9, 
attain a maximum value near p = 1.04 and decrease there- 
after. 
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TABLE 3.1 


Seme Approximate Solutions 

P = 0.05, R = 0.2, = 0.02, ^2 = 0*0^? 

= 0.5 


Frequency 

P 

X Number of | 
i harmonic x 

I solutions I 

ij 

e - amplitude | x - 
B 0 

amplitude 

A 

X stability 
i characte* 
^ ristic 

0,9^ 

1 

1.48 

0,23 

stable 


1 

1.62 

0.30 

stable 

0.96 

2 

2.33 

1.00 

unstable 


3 

2.53 

1.59 

unstable 


1 

1.79 

0.40 

unstable 

0.97 

2 

2.06 

0.61 

unstable 


3 

2.50 

1.47 

unstable 

0.98 

1 

2.45 

1.31 

unstable 
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Lastly Fig. 3,5 shows that for all -/alues of P, the 
solutions become unstable v/hen amplitude of Xj A, reaches 
a value near 0.35 making it the critical value of A for 
the combination of parameters selected. 

For representative points lying in region II , Fig. 
3.3} the amplitudes of the corresponding unstable solutions 
are not large as they were for region III. In fact, it 
can be seen from Figs. 3.^ and 3*5 that as region II is 
entered while decreasing p, the amplitudes, A and B, 
decrease. Consequently there is not much difference between 
the results obtained by successive approximations. Only the 
boundary of region II is sli^tly dependent on the order of 
approximations . 

Finally typical plots of phase angles, and 2 
with P =0.07 are shown in Fig. 3.6. From this figure it is 
observed that as region III is approached with increasing p, 
the phase angles also change very sharply. Similar trends 
are obtained for other values of P as well. 

It may be noted here that for q = » 

the linear natural frequency of the pendulum taken as half 
of that of the primary system) the region III starts at 
p 1.05 with P = 0.05, R = 0.2, = 0.02 and 
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Another jnteresting point is that the transition 
from stable harmonic solution to unstable harmonic solution 
of region III is associated with sharp increase in the 
amplitudes 5 viiereas, the transition to region II is smooth* 

3.5.3 Numerically Integrated Results 

In this section the numerically integrated results 
are presented. These results are compared with the approxi- 
mate results discussed in section 3.5.2. First the boun- 
daries of regions of Fig. 3.3 are considered. Then the 
results for regions I, II and III are presented sequentially. 

3 . 5 . 3 .1 Boundaries of Regions of Stable and 
Unstable Harmonic Solution 

The boundaries of regions of stable and unstable 
harmonic solutions are obtained by numerically integrating 
(3.1) and checking whether the solution reaches a harmonic 
(or near harmonic) steady state or not. For different 
regions, the boundary points by integrations are found by 
increasing or decreasing p in steps. This is done because 
near these points the responses are strongly dependent 
on initial conditions. Boundaries of these stable and 
unstable regions are sho^ci in Fig. 3.7. Sane results 
obtained by the approximate method, already ^own in Fig, 3.3, 
are also included. The boundary points frem integrated 
results are indicated by cross marks. This figure reveals 



BOUND/iRIBS OF STABLE AMD TUI STABLE HAOfCHZC 
SOLOTICKS Bt APPHQXIMATB METHOD AND 
iSTBGRATiraiS 








that for O.Oif < p < 0.09 the integrated results yield 
a larger span of the left hand side overhang than that 
obtained by approxims-te solutions. For integrated results 

at P = 0,09 region II widens sharply towards left and 

< 

touches the boundary of the left overhang. After this the 
left overhang span almost disappears, i.e,, with P > 0,10 
the left boundary of region II almost coincides with the 
boundary of primairy instability. It is evident from Fig, 

3.7 that, except for the left overhang, the results obtai'ned 
hj integration and approximate method match very ^;ell, 

3.5* 3» 2 Numerical Results for Region I 

The amplitudes of 6 and x are shown in Figs. 3*8 
and 3 .9} respectively, where both the integrated and the 
approximate results are included. For the integrated results, 
when the responses are seen to be periodic (.with x and 6 
having periods of 2% and respectively), the maximum 
value in one period of oscillation is taken as the represen- 
tation of the amplitude. For the cases, vhexe tiie amplitudes 
obtained by the approximate method and integration match, 
the integrated responses arc, of course, seen to ce 
completely harmonic. For representative points in region I 
of Fig. 3 , 7 , the integrated responses reach harmonic or 
near harmonic steady state. Hiese steady state amplitudes 
are joined by firm lines in Figs. 3*8 and 3.9* dotted 

lines in these figures are for ccaabinations of P and p 
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•FIG, 3*8 


P 

liWliTSWM OF PHfDULiOM OBTAINED BY APPROXI- 
MATE MEEHOD AND INTBGRATICK 


i 







FTG <? AMPLITUDE OF PRIMAHI (BTAINED BT 
FIG. 3.9 ^JrqximATB method AND INTEGRATICK 
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lying in region II of Fig, 3,7, For these values the 
integrated responses do not have fundamental period of x 
and e as 2 % and 4 71 , respectively, and the amplitudes of 
responses do- . not have any meaning in ordinary terms. In 
fact, the dotted lines in Figs. 3.8 and 3*9 are nothing 
hut the extrapolations of the respective firm lines. Figure 
3*8 shows that the amplitudes of 9 ohtainad hy the approxi- 
mate method and numerical integration are in very close 
agreement. At lower values of p, however, these two methods 
give different values for the amplitude of x (Fig, 3.9)* 

The reason for this discrepancy is explained in the follow- 
ing paragraph. 

Figures 3.10 and 3.11 sha-; the integrated steady 
state response time histories of D and x, respectively. 
These responses are for P = 0,07 at different values of p. 
Figure 3,10 shovrs that the 0 - response is harmonic for all 
values of p. Thus the assumed foim, given hy (2.l6), is 
correct and so the results obtained hy the approximate 
method are very accurate. Figure 3.11 shows that in the x - 
response, the contributions of the higher haimonics increase 
with decreasing p. This explains the error in approximate 
solutions where only the first harmonic was consiaered. 

More the contribution of the higher harmonics, greater is 
the discrepancy between the approximate and the integrated 


results in Fig* 3.9. 



21 ? 



FIG. 3-10 
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V/e now discuss some results for P= 0 .l 5 , In 
region I, both the approximate method and the integration 
give stable harmonic or periodic solutions only for two 
very narrow bands of p — one around p = 0,84 and the 
other around p = 0,90, The typical solutions at p = 0 , 84 - 
are: X = 0 , 34 - cos ( t - 0,22), Q = 1,96 cos ( ) by 

approximate method and x = O.32 cos ( t - 0.20), 
e = 1»95 cos ( 2 ^* 92 ) by integration. At p = 0,90 the 

approximate method yields a stable harmonic solution: 

X = 3.81 cos ( t - 1 , 88 ), e = 3.31 cos ( I + 0.61). The 
integration gives a strongly stable steady state as 
X = 3,85 cos ( t - 1 . 4 - 4 -), e = II + 3 . 64 - cos (| - O.76). 

The amplitudes of x and 0 are seen to be very high with 
that of 0 being > % , 

3.5*3 *3 Numerical Results for Region II 

For points in region II of Fig, 3 . 7 } the integrated 
responses of x and 9 are found not to have the fundamental 
periods as 2 % and 4 -ii , respectively. The responses may 
or may not be periodic depending on the initial conditions 
and the values of P and p. Moreover, again depending cn 
initial conditions, more than one periodic steady state may 
exist. Some idea about the existence of periodic or near- 
periodic responses can be obtained from the eigenvalues of 
(2.32) used in the approximate method. This is explained 


below with sane examples 
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Table 3*2 gives the eigenvalues of (2.32) for three 
different cases. In each case, the value of A, B etc. of 
(2.15) and. (2.16) are also given In this table. The approxi- 
mate solutions of steady states are used to generate the 
following set of initial conditions by taking t = 0 in 

( 2 % 15 ) and ( 2 .l 6 ); 

X ( 0 ) = A cos 

X* (0) = -A sin $ . i 0*12) 

' > 

0 ( 0 ) = Be os §2 

e'( 0 ) = - I sin §2 

Figures 3.12 to 3.l4 show the integrated steady state 
responses of 0‘for the cases cited in Table 3.2. The res- 
ponses are presented for 0 ’ instead of 0 because in some 
cases, even in the steady states, the pendulum is found to 
make complete revolutions, periodically or non -periodically. 

m case (i) (Table 3*2), P and p ccmblnation is in 
region I but is very close to region II. The maximum 
eigenvalue is just less than unity. Therefore, the steady 
state given by the approximate solution is expected to be 
weakly stable. Figure 3.12 Xb) shows that with initial 
conditions given by ( 3 . 12 ), 0 * - response is haimonic with 

period ifTc . However, with slightly different initial con- 
diticxis, the integration carries the solution to a different 
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TABLE 3.2 

Eigenvalues of (2.32) for Three Different Cases Near 
Region II 


Case 

1 Parade ter 
A comhina- 
1 tion 

1 

1 Amplitudes 

A and phase 

5 angles of 

0 (2.15) and 

5 (2.16) 

! Eigenvalues of 

(2.32) ^ 


y 

t ^ 

S 1 

5 

i) 

S 

o>»xr>«»<o 

C\J 

¥ 

3 s 
j 

‘if 



A = 0.098 

-0,51 + 


-0.51 - 

-0.28 + 

-0.28 - 


P = 0.06 

B = 

0.85 i 


0.85 i 

0,46 i 

0.46 i 

(i) 









p = 0.88 

<{)^ — 0«63 

modulus 

== 


modifLus : 




<l>2 = 0.01 

0.99 



0.54 




A = 0.12 

-0.70 + 


-0.70 - 

-0.32 + 

-0.32 - 




0.78 i 


0.78 i 

0,40 i 

0.40 i 


P = 0.07 

B =1.35 






(ii) 









p = 0.88 

= •"0.51 

modulus 

z: 


modulus : 

= 



= 0 

1.05 



0.51 




A = 0.12 

-0.93 + 


-0.93 - 

-0.33 + 

-0.33 - 


P = 0,08 

B = 1 .46 

0.75 i 


0.75 i 

0.28 i 

0.28 i 

(iii) 









p = 0.86 

= -0.49 

modulus 

z= 


modulus 




<l>, = 0.02 

1.19 



0.43 
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steady state as shovm in Fig, 3*12 (a). Hie response in 
this case is scon to be amplitude modulated with a time 
period of 127c . 

For cases (ii) and (iii) (Table 3*2) the representa- 
tive points lie in region II and the maximum eigenvalues are 
slightly greater than unity. The corresponding steady statesj 
obtained by the approximate solutions, are therefore unstable. 
For numerical integration also no initial condition is found 
which takes the solution to a harmonic steady state. As the 
mazimum eigenvalues in these cases are complex and near 
unity it is expected that the growth of the trivial solution 
of (2.27) and (2.28) will be slow and oscillatory. There- 
fore, with the initial conditions- (3.12), the numerical 
integration is expected to give a steady state which will be 
close to the harmonic steady'' state given by the approximate 
solutions. This is confirmed in Pigs. 3.13 and 3.1^j which 
show the Integrated steady state responses of e* with two 
different sets of initial conditions. The firm lines in 
these figurec correspond to initial conditions (3.12). 
Comparing Figs. 3*13 and 3.1^ it is seen that the responses 
become more irregular as the maximum eigenvalue increases. 

The dotted line in Fig. 3*13 shows another steady 
state e’ - response for case (ii), obtained with a different 
initial condition. Though the trends are similar, this is 
seen to differ slightly In magnitude from that represented 



87 


tiy the iirm line. It may to noted that both the responses 
are periodic v/ith time period 44 ic . 

Similarly the dotted line in Fig. 3 *14 shows another 
steady state obtained with an initial condition which is 
different from that used for firm line. The fiirni line shows 
a period of 20 u , whereas, the dotted line response has a 
period 2 jx , Moreover, in the latter case, the mean value 
of 0‘is nonzero. This implies that during the steady state 
the pendulum keeps, on rotating in the same direction, the 
time for each rotatioi being 2 % • 

The initial conditions, mentioned above, yield 
periodic steady states in each ease considered in Table 3»2. 
However, if initial conditions are chosen so as to be 
sufiiciently far away from these periodic steady state 
tra;iectories, the integrated response becomes irregular and 
nonperiodic • 

V/ith further increase in P than mentioned in Table 
3»2 (with the representative point still lying in region II) 
the maximum eigenvalue also increases and the integrated 
responses become nonperiodic with any choice of initial 
conditions . 

All the above mentioned results are verified with 
a smaller step size of integration ( = ^ ). 
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3*5.3*^ N-ujiierical ResiiLts for Region III 

For points in region III of Fig. 3 •7? vfhile checking 
the stSihility of the approximate harmonic solutionj the 
maximim eigenvalue of (2.32) is found to be real and much 
greater than unity. This implies that no harmonic or near 
harmonic steady state can be obtained in the nei^bourhood 
of the solutions obtained by the approximate method. In 
such cases the integrated responses are seen to be non- 
periodic. In other words, even after a long time, when the 
effects of the transients are expected to die down, the 
responses do not show any periodicity and are quite irre- 
gular. Moreover, the results do not converge by decreasing 
the stop size of integration, 

' For points in region I, but very close to region III, 
harmonic steady states are obtained only wilii tiie following 
scheme: The integration is started with points well inside 
region I and lAhen the steady state is reached, the frequency, 
p, is increased in small steps. The harmonic steady 
states, thus obtained, arc very weakly stable. As an 
example, we take P = 0,05 p = 1*0, when the approxi- 

mate method gives an unstable harmonic solution with the 
maximum eigenvalue being real and much greater than unity. 
Following tie scheme merticned above, a harmonic steady 
state, with x = 0,3 cos ( t - 0,6) and 6 = 1*52 cos 

(■| - 0.2) is found by integration. Hcn.rever, this steady 



89 

state is very weakly stable so that even a disturbance of 
0,1 to the coordinate e , takes the solution away from this 
steady state and the responses become nonperiodic. For all 
other initial conditions; the responses are completely 
nonperiodic and the results do not converge oven with step 
size equal to . Figure 3-15 shows typical time his- 
tories of e‘ for txiTO different step sizes ( 9(0) = 0,05, 

e* (0) = X (0) = x’ (0) = 0). 

For points well within region III no harmonic or 
periodic steady state can be obtained. For one such point, 
with P = 0.05 and p = 1.05, e’ - response is sho^m in Fig. 
3.16, Figure 3,17 shows x - responses for P = 0.05, 
p = 1,0 and 1.05. Since x - coordinate is directly driven 
by the external harmonic excitation, it is expected that 
X - response will have more regularity as comparea to that 
for 0 and this is revealed by Fig. 3 *17 as well. 

The results, presented in Figs. 3.15 to 3-17, as such 
do not convey any meaning since these results have not con- 
verged with decreasing step sizes of integration. Numerical 
calculations in double precisions, cmpl eying better scheme 
of integration by predictor-corrector method, also fail to 
converge. Moreover, a change even in the eighth decimal 
place, in the initial conditions changes the response com- 
pletely. These differences in the results for slightly 
different initial conditions occur because it is found that 
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scmewhere In the tine history, 9 ’has a veiy steep gradient. 
This makes the subsequent ccmputations very sensitive to the 
inlierent inaccuracies of numerical integrations. In compa- 
rison! for region I, where the harmonic steac3y states are 
strcngly stable, even a coarse stop size of integration with 
any initial condition carries the solution to the sane har- 
monic steady state. 

To derive any meaning out of tiie results obtained for 
region III, the following argument is put forvra.rd. In the 
actual physical system, if experiments are performed then at 
every stage there will be some unaccountable external dis- 
turbances. It is also not possible to control the parameters 
to exact values. The numerical integrations can be thought 
of as a computer simulated experiment. Here, in fact, the 
control on the parameter values is much better than in the 
actual physical experiment. The external disturbances in 
the physical system may be represented by the numerical 
round'^fe and the inherent errors in the numerical schemes 
used. 

In the light of the above arguments, a strongly stable 
solution implies consistent and reproducible results from a 
physical experiment. Likewise, an unstable or a weakly 
stable solution suggests that the correspending experimental 
result can not be reproduced. This is because of the inhe- 
rent disturbances which, in the latter case, change the 



response drastically. These disturbances can only bo consi- 
dered as random. So it is iiapossible to predict the exact 
value of the response at a future instant of time. The 
response time histories can then be treated as sample records 
of a random process, 

"In practical terms, the decision as to whether or 
not physical data aio deterministic or random is usually based 
upon the ability to reproduce the data by controlled experi- 
ments [4-8]". For parameters lying well inside region I an 
experiment, producing identical results, can be repeated many 
times (within the limits of experimental error) and the res- 
ponses, thus, are deterministic. For region III the experi- 
mont can not be designed which will product identical results 
when the experiment is repeated and then the data must be 
considered random in nature. As discussed above, the inte- 
grated responses can be considered to represent seme physical 
experiment. Therefore, the only meaningful description of 
'these time histories can be attempted through the statistical 
analysis, which is presented in Chapter 4-, 

3. 5. 3*5 Numerical Eosults for Sufficiently High Excitation 

It is clear from Fig, 3.2 that if P is large then 
for certain frequencies the representative point lies in 
the secondary region of instability . For example, with 
P = 0,20 and p = 0,90 the numerical results give a periodic 
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steady state with e oscillating with period 2i; . e -* 
oscillations with a period of 2n are usually seen to he 
associated v/ith a shift in nean values of 0 from zero* 

This shift of the dynamic mean position from the static 
equilihrium position is well known in nonlinear systems and 
is often called drift or steady - streaming [9]* 

With sufficiently large P while increasing p the 
secondaiy and the primary regions are alternately trajrer^ 
sod in Fig. 3 •2. As an example, with P =1,0 the follow- 
ing behaviours of the steady states are observed while 
increasing frequency: (i) For p <0,30 the zero solution 
of e is stable, (ii) At p = 0.4-0 the representative 
point is in the secondary region. 9 - oscillation has a 
period of 2% and a mean value of 0,23. (iii) At p = 0.55 
the zero solution of e is a stable solution, (iv) With 
p = 0.70 the point is in the primary unstable region and e 
oscillates with a period of 4m . e and e’ responses are 
shoxm in Fig, 3*18 (a)* (v) Increase in p to 0,80 takes 
the point to secondary -region again and © oscillates with a 
period of 2m . Figure 3*18 (b) shows e and e’ responses 
and the mean position of 6 - oscillation is near g • 

(vi) At p =0,91 e oscillation is seen to have a period of 
2 Tc and therefore the point, perhaps, lies in the foiirth 
region of instability, the boundaries for ■v^ich are not 
computed, (vii) At p = 1.0 the point is in the secondary 
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unstable region and 9 oscillates with a period of 2 % • 

(viii) For p greater than 1,20 the pendulum is in the 
primary region of instability. At p = 1,20 and lAO the 

pendulum is found to make continuous rotations with a period 
2 % j whereas, at p = 1 ,80 it oscillates with a period of 
^Ti: • (ix) Finally at p = 1,90 tho static equilibrixim 
position of the pendulum again becomes stable. 

Above are the cases at certain frequencies where x 
and 0 oscillate periodically. In between those frequencies 
there are regions where the integrated responses are non- 
periodic . 


CHAPTER if 


STATISTICAL ANALYSIS OP THE RESPONSE TBIE HISTORIES 
if.1 INTRODUCTION 

In this chapter the nonperiodic responses j obtained 
by numerical integration, are analysed statistically. Some 
of these responses were sho>j:a in Pigs. 3.1 5 to 3.17. All 
the relevant standard definitions and formulae used in this 
chapter are given in appendices A-1 and 2, respectively , 

The oustification of doing the statistical analysis was 
given in section 3.5»3.^. It was also mentioned there that 
the responses for the nonperiodic cases are dependent on 
the initial conditions and the step sizes of integration. 

It is assumed now that with other parameters rmaining same, 
the integrated response record, for each combination of the 
step size and the initial conditions, is a sample record 
from the same random process. Thus by changing these two 
parameters different sample records are obtained. 

It is further assumed that these sample records 
are the realizations of a weakly ergodic random process. 

In other words the sample records are considered uo be 
stationaiy. The assumption of ergodicity renders the faci- 
lity of taking time averages over a specific saii 5 )le record 
instead of ensemble averages. 
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First the sample records are tested for the 
hypothesis of stationarity • This test is done throu^ the 
run tests on the mean sq.uared values calculated over diffe- 
rent intervals fraa the same sample record. Once this 
hypothesis is established satisfactorily, different statisti- 
cal descriptions like autocorrelation functions and power 
spectral density functions (PSD) are obtained. The near 
constancy of these quantities over different records 
(obtained with different combinations of the step size and 
the initial conditions) justifies ttie assumption of treating 
each response hist-ory as a sample record. 

The statistical analysis is done for the sample 
records of 5E and 6 * . ©'coordinate is selected instead of 
0 » This is done in order to avoid confusion as in the 
nonperiodic cases, the pendulum is found to make sane can- 
plete revolutions in either direction without any fixed 
pattern. 

The statistical quantities ^own to depend on 
the parameters P and p. in tnis work no attempt is made 
to correlate the statistical quantities with P and p* 

^,2 PAHAMETER ESTBIATORS FOR THE SA14PLE HB30BD 

For a stationary sample record {u^} j taken at 
N points, with the time intei^al between the successive 
points being h^, the estimators for the basic descriptive 
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properties (defined in appendix A-2) are given telov;: 
(i) Mean value estimator ( {1 ) 


u 


1 ^ 

N 2 

i=1 


(4-.1) 


For ease of calculations, the data points at this stage are 
transformed as 


u„ - P-. 




; n = 1, 2, 

'n n u ’ 7 j 7 

The following estimators are given for tlae transformed 
record { v^} . 

Cii) Variance estimator (cf^ ^ 


(if.2) 


^ 2 

'^v 


1_ 

W - 1 


N 

2 

i=1 


(^•3) 


Equation (4.3) gives the unbiased estimator for . 

Since the data are transformed, (4.3) also gives the mean 

square value. 

A 

(iii) Autocorrelation function estimator ( 

■5 _ __J ^2^ V V „ 5 3? = 0,1,2 , ...,s, 

R^ = K ^ r ^ n n+r ’ 

i-1 (4.4) 

where s is the maximum lag number. R^ is the unbiased 
estimate of the autocorrelation function at lag r. 

The normalized autocorrelation function R^ is 

A ^ 

obtained by dividing R^. by R^ , i«e.j 



1C1 


R 

R = ; r = 0, 1,2, ,.,,s 

^ R 

o 


(4.5) 


The normalized autocorrelation fimction for x and 
0* are represented by R_ and R j , respectively (dropping 

X ® 

the subscript r). Variations of these normalised auto- 
correlation functions with the lag' number, r, are shown in 
the autocorrelograms. 

(iv) Power spectral density function estimator ( ' 


PSD estimator is obtained via the normalized 
autocorrelation function estimate in two steps. First the 
raw estimate is obtained as 


s-1 

G, = G (f ) = 2 h, [ 1 + 2 S R cos ( ) 

r=1 




V 


+ R COS J 5 

s 


k-rt ] 


(4.6) 


where G, is calculated at m + 1 discrete frequencies 




k= 0,1,2, *•*, 


(4.7) 


with f as 
by 


the cutoff frequency (Nyquist frequency) given 


(4.8) 


2 h4 


The smoothened PSD, is then calculated 


as [ 46J 
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A 



0*5^ + 0.4-6 G, 

o -] 


0.23 + 0.54 Oj^ + 0.23 i 

0.46 G^_^ + 0.54 G^ 


^ S*^ *1 

(if.9) 


Dropping the suhscript k, the smoothened PSD's for x and 9 * 
ai e represented hy G_ and Gg* j respectively. 

^.3 TEST FOR STATIQNARITI [48] 

The concept of stationarity, defined in appendix 
A~1, relates to the ensemble averaged properties of a random 
process. However, in this chapter the individual sample 
record, assumed to be from an crgodic process, is referred 
to as being stationary. By this it is meant that, except 
for the variations due to normal statistical sampling, the 
properties computed over short time intervals do not vaiy 
significantly from one interval to the next. We assume that 
the time trends in the mean square value will reveal the 
nonstationarity of the sample record. Moreover, the sample 
record is considered to be long canpared to that demanded by 
the lowest frequency component in the data. It is further 
assumed that the sample records frcm a nonstationary process 
will also be nonstationary. With these assimiptions tiie 
stationarity of the random proces| is tested by considering 
a single record and therefrom compisting the variances over 
short time intervaJ.s and'peri&SJming "the run test given in 


appendix A- 3* 
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h.h RESULTS MD DISCUSSIONS 

All the result^ presented in this section are for 
the f olio-'-' ig values of the parameters: R = 0,2, ^.j = 0,02, 

^2 ~ 0.2 ~ ^ 0 = 0-]= 0,5* The parameters varied 

are P and p. Since the nonperiodic responses depend on the 
initial conditions and the step size of integration, two more 
parameters are taken into consideration, namely, initial 
c-ondition of 0 , i.e., 0(0) and the step size of integration 
h. The initial conditions for other coordinates are taken 
as zero. 

It 

The data points are taken at interval ~ 10 * 

A further decrease in the internal size does not change the 
results. Moreover, the cutoff frequency, f^ , with h.^. = 
is i- ( Eqn, ^.8) which is quite large as ccmpared to the 
highest frequency components of the sample records as shown 
later. If H data points are taken the total record length 


i+.i+.l Mean Values and Variances 

Mean values of x and ? , for all the sample 
records considered, 'are found tc he almost zero. The follow 
ing results are for the transfoimed data. The variances 
are oaloiiLated with different nnnters of data points (l.e., 

different record lengths) ty using (4.9). lables 4.1 (a) 

— ^ 2 ^2 

and .^+.1 (h) show the variances of x and 6 , <y£ and 0^', 
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TABLE 4.1 (a) 

Variances with Different Ilunher of Data Points 


if 

o 

• 

P 

= 1.0, 

0 (0) = 

0*05j h 

Nmbe of data 


Record 

5 

Variances 

points 

N 

5 

0 

lengtli 

5 A 2 

S 

8 * ® 

! «e' 

200 


20^ 

0.077 

0.19 

400 


4011 

0.064 

0.33 

600 


60 % 

0,065 

0.31 

800 


80 It 

0,084 

0.34 

1000 


100 It 

C.079 

0.35 

1200 


120 71 

0.078 

0.34 

1600 


160 It 

0.079 

0.35 

2000 


200 7C 

0.079 

0.36 

3000 


300 t 

0.076 

0.35 

4000 


400 % 

0.074 

0.33 

8000 


800 % 

0.077 

0.34 


10 ? 


TABLE 5+.1 (b) 


Variances witii Different Humber of Data Points 

p = 0 * 05 , p = 1.05 e(o) = 0.05, h = 1^:5 


H'miber of data 
points 

H 

5 

il 

Record 

length 

r ~ 

2 

2 

2 

Variances 

2 * 2 

2 

2 

200 


20 11 


0.092 

0.29 

ifOO 


5f0 % 


0.067 

0.21 

600 


60 % 


0,065 

0.21 

800 


0 

00 


0.056 

0.18 

1000 


100 n 


0.053 

0.20 

1200 


120 IX 


0.050 

0.20 

1600 


I6O ix 


0,048 

0.19 

2000 


200 Tx 


0.054 

0.19 

3000 


300 7X 


0.056 

0.20 

5+000 


5+00 IX 


0.053 

0.20 

8000 


800 tx 


0.055 

0.20 


PGspectively , fop tyo cases uith increasing value of N, It 
can be seen that the variances do not change significantly 
for N > 1000, 

^»4-.2 Run Tests for Checking Stationarity 

Following the procedure outlined in appendix 
each sample record is divided into 20 equal time intervals 
with each interval having 200 data points (i.e.^time interval 
= 207r )* For each interval the variances are calculated and 

then these are aligned in time sequence to perform the run 
test • 

Figure ^,1 shows the values of the variances versus 

the interval number for two combinations of parameters. In 

each case the median value is shown by the dotted horizontal 

line. V^ith respect to this median value, tlie number of runs, 

r , is counted. It is seen that r„ lies between 6 to 15. 
n ’ n _ 

Thus, from the table in appendix A-4-, the hypothesis of 
stationarity of the sample record is found to be acceptable 
with 0,05 level of significance. 

Similar run tests are perfoimed for other caabina- 
tions of parameters. Seme of these values are given in 
Table ^.2, where only the median values and the number of 
runs are indicated. Table 4-*3 shows similar results with 
longer time intervals, each interval now having 1000 data 
points (i.e.jtime interval = 100 u ). This Interval size 
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TABLE h,2 


Median Values of Variances and Number of Buns 


Number o' 
interval s = 20 


Number of data points 
in each interval = 200 


{ { } { 

loi I M 

5 ^ ^ ^ - 


Median Values of 0 Number of runs of, 
I~2 5 * 2 X 

ff - ^ aj 5 s {J„r 

^ 1 ® 8 ^ il ^ 


1 

0.0? 

1.0 

0.0? 

n 

0.07if 

0.28 

12 

9 

2 

0.0? 

1.0 

0.0? 

% 

m 

0.082? 

0.30 

7 

12 

3 

0,0? 

1.0 

0.10 

n 

5o 

0.0?79 

0.33 

11 

11 

if 

0.0? 

1.0? 

0,0? 

u 

5o 

O.Oif? 

0,16? 

11 

11 

5 

0.07 

1.0 

0.0? 

% 

^0 

0.09?9 

0.if02 

11 

10 

6 

0.07 

1.0? 

0,0? 

n 

5o 

0.071 

0.223 

11 

7 

7 

0.07 

1 .08 

0.0? 

TC 

5o 

0.063 

0.21 

12 

11 
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SI. 5 $ 

No. 5 


I P 


8 Mndian values of j Number of runs of 

^ H2 r~*2 f 71 J * 2 

U O <j^ cj- I 0^, 

5 ^ } i) ii ^ 


1 

0.05 

1.0 

0.05 

0.076 

0.3^9 

14 

12 

2 

0.05 

1.0 

0,10 

0.075 

0.342 

12 

10 

3 

0.05 

1.05 

0,05 

0.055 

0.205 

13 

11 


0.07 

1.0 

0.05 

0.095 

0.45 

13 

■10 

5 

0.07 

1-05 

0.05 

0.0817 

0.37? 

11 

9 

6 

0.07 

1.08 

0.05 

0.062 

0.228 

1? 

17 
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is tested because it was seen in Table 4-.1 that the variances 
do not change nuch for N > 1000. It is obvious from these 
tuo tables that only, in the last case of Table tile 

number of runs for ^ barely exceeds 15* Several other 
sample records were also tested. For all the cases, except 
the one mentioned above, the hypothesis of stationarity of 
the sample records is found to be acceptable with O.05 level 
of significance. Thus with the hypothesis of stationarily 
having been established, further statistical analysis is 
undertaken in the following sections. 

^.4.3 Autocorrelation Functions and Po\/er Spectral I 

Density Functions 

Normalised autocorrelation fimctions and smoothened 
power spectral density functions are calculated using (4,5) ! 

and (4,9), respectively. VJith time interval h^ = -riQ the 
cutoff frequency ! 

f = 5 (4.10) ! 

c % 

Fran (4,7) the PSD functions are calculated at the following ! 

(s + 1) discrete angular frequencies: i 

I 

; k = 0,1,2, s (4.11) I 

K. S I 

I 

Figures 4*2 and 4,3 show the typical autocorrelo- i 

I 

grams for x and '9, respectively, with P = 0*05, p = 1 *0, ! 

0(0) = 0.05 and h = * These are obtained by taking ■ 

i 

1 
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^-000 data points (l.e.,N = 4000), The maxinian lag number, 
s, Is taken to be 10^ of N, i.e.,s = 400, Almost similar 
curves are obtained for N = 1000, 2000 arid 8000. For the 
same ccmbination of parameters and with different values of 
l‘I and s the PSD functions for ic and e’ are shown in Figs. 

4,4 and 4,5. It is seen from these figures that with 
s = 0,1 N, different values of H (i.e., different record 
lengths) give the same trend for PSD's. Some results are 
also shown with s = 0,2 H and s = 0,05 N. It may be noted 
that so far as the PSD functions are concerned, sufficiently 
accurate results are obtained with II = 1000 and s = 100, 
Further results, unless mentioned otherwise, are obtained 
■with these values of N and s. 

For the sane values of the parameters used in 
Figs, 4,2 and 4,3, with the step size only changed to ^ , 
the autocorrelograms for x and 6* are shom in Figs, 4,6 
and 4,7 j rospectively . The trends of R_ and P-gtin these 
figures are seen to be similar to those of Figs. 4.2 and 4.3, 
For the same values of P and p and for five different 

A A 

ccmblnations of 0 (0) and h, G and G , are shovin in Figs. 

X e 

4,8 and 4,9, respectively. It may be concluded from these 
figures that -the different cembinatiens of 6 (0) and h 
(i.e., the -independent realisations) yield the same trend d 
PSD’s implying that all the realizations are from the same ‘ 
random process. 
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h»h-,h Response Statistics 

Now having established that for a given conbination 
of P and p, any value, of 0(0) and a proper value of h 
yield a correctly representative stationary saii'ple record, 
seme response statistics for different conbinations of P and 
p are discussed below. 

(i) P = 0.05, p = 1.0 

With this parameter conbination the representative 
point lies near the boundary of regions I and III in Fig. 
3,3* It may be recalled tliat near tills point there Vas a 
sharp increase in the amplitudes (Figs. 3.^ 3.5)* 

Figure ^.3 shows that N0» , almost tends to zero 
as the lag number, r, increases. Hiis suggests that e’ 
behaves like a randon variable. Frcra the nature of Rgi it 
may be further observed that the sample record of behaves 
like that of a wide band random process. On the other hand, 
Fig. ^.2 shows the presence of a strong periodic trend in 
the record of x. However, as R_ falls below 0,5 for 

X 

!> > 20, the data can not be considered as deterministic. 
Piather it can be taken as a randan noise superimposed on a 
periodic signal. The frequency contents of the responses 
are better seen in Figs, 4- ,4, 4.5, 4.8 and 4.9. 

Figure 4.8 shows a large value of PSD for x at 
the non dimensional angular frequency m = 1.0 (i.e. at the 
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forcing frequency)# This is further confirnod in Fig. 
where, hy increasing the record length, the peak at u = 1 ,C 
increases almost proportionately. Ihe presence of a strongly 
periodic component in x - response at the forcing frequency 
is also expected hecausc the primary mass is directly driven 
by the external excitation. The time history records of x 
in Fig. 3*17 also show a predominant periodic component. 

Figure 4-. 9 reveals low frequency contents for 
6 - response. Ho peak is observed at cu =. 2 » ^*7 

shows that PSD characteristics remain unchanged with increase 
in record lengths. The corresponding time history was shown 
in Fig. 3«15» 

(ii) P = 0.05, p = 1.05 

For this conbination the representative point is 
in region III (Fig, 3 , 3 ) and lies near the boundarj'- of primary 
instability . 

PSD of X - response is found to be similar to that 
of Fig. 4-. 8 and is not shown. Autccorreloeram of © ‘ - response 
is shown in Fig, 4-, 10, It is seen that the data now have 
a periodic trend. PSD of s' - response (Fig, 4-, 11), besides 
having a low frequency content, also has a peak at a = 2 • 

By increasing the record lengths this peak is seen to 
increase almost proportionately, whereas, the low frequency 
contents remain almost same. Tlius, the sample record of q* 
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— -1 

has a periodic coiiponent at oj = ■^. The corresponding time 
histoiy v/as shown in Fig. 3*16 which looks regular caapared 
to that in Fig* 3»15* 

Ciii) P = 0 * 07 , p =1*0, 1.05 and 1*08 

For p = 1.0 and 1*5 the representative point is 
well inside the regien III in Fig. 3*3* At p = 1.08 the 
point is close to the boundaiy of primary unstable region. 

In each case the sample record is obtained with 
h = and 0(0) = 0.05. PSD*s for x and 9* are shov/n 
in Figs. 4.12 and 4.13, respectively* As in the previous 
cases, PSD for x shows a peak at (1 = 1. For p = 1.0 and 
1 . 05 , PSD of 8* has only low frequency contents without any 
predominant peak* As p is increased to 1.08, a peak is 

— -J 

observed at &) = -^ , 

For all cases discussed in this section, the median 
values of the variances are given in Table 4.3* For non- 
periodic responses, the variances together with the PSD*s 
should be ijsed to describe the response statistics# 
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CHAPTER 5 


©IPERIl'IENTAL EIVESTIGATION OF IIOl.’PERIODIC EESPaiSES 

5.1 BTTRODUCTI® 

In this chapter sene experimental studios are 
presented for a ha^rmonicaily excited t\/o degree- of -freedai 
autoparametric system. The 33’"stem, considered in this 
chapter, is similar to the ones studied in the previous 
chapters except that the system now is base-excited instead 
of being force- excited. Though the eauaticiis of motion for 
these two systems are slightly different, the basic charac- 
teristics remain same. 

It is seen in chapter 3 thar for certain ccinbina- 
tions of parameters 1iie numerical integrations yield non- 
periodic results. In such cases the nonconverging results 
are later analysed in chapter 4, based on the stipulation 
tliat these still represent some physical phenamena where the 
corresponding physical experiments will also yield non- 
reproducible results. In this chapter some experimental 
investigations are undertaken to justify the above mentioned 
hypothesis qualitatively. The parameters in the experiment 
are adjusted so as to cause instability of the harmonic 
steady state motion giving rise to random responses. The 
statistical analysis of the experimental results is not 
undertaken . 
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5.2 MJIEIEMATICAL MODEL AI^D EQUATIONS OF MOTION 

Tho aathematical model of the experimental system 
is shown in Fig, 5.1. The primary mass, M, is connected 
to the base ^irith a linear spring of stiffness and a 

■viscous dashpot with damping cons'bant c’^. A cempound 
pendulum of mass m is hinged to -the primary mass. The 
effective length of the pendulum is 1 and its moment of 
inertia about the controidal axis is I. The pendulum ro'ba- 
tion is ass'umed to be damped by a viscous torque for which 
the damping constant is 02 * The base displacement is repre- 
sented by y - coordinate and the absolute displacement of 
the primary mass is sho-wn by x; e is the angular displace- 
ment of the penduliM. 

The eq’uations of motion of the system are 
(M + m) X + c^ ( X - 3 O + (x - y) 

- ml ( i sin e + w cose ) = 0 (5.1 ) 

and (ml^ + l)e + C 2 e + nil (g - x) sine = 0 (5.2) 

The base is given a harmonic excitation 

y = Y cos “t , (5*3) 

■vdiere Y is the amplitude and « is the angular frequency 
of excitation. Using (5.3)j (5*1) a^id (5*2) can be written 
as 



r, 

r 

’x 


.FIG, 5 * 



MATHEMATICAL MODEL OF IHE EXPSE! 


I'ljffiSTAL 
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x+2 uj^x + to^x - ( e sin e + cose ) = 


4 


and 


Y / {D^j + (2 03^ 03 ) COS [ <at + tan*" ( 

V 


1 2 a 

-1 ^ L,) ] 


a 


1 


0+2^2 “ 2 ® + C“2*'^f ° » 

respectively, whore 


(5A) 


( 5 . 5 ) 


05^ 


/ 

= — *— is the natural frequency of the ^ 


locked mass system, 


is the inertia factor, 

ml'^ + I 


(iJ, 


= ^ ^ ^ is natural frequency of the 

ccanpo-und pendulum, 


^5*6) 


R = 


E 

11 


is the mass ratio, 


- 

C. = 7 ^ TP T is the damping ratio of the 

1 2 vli + ml r 0 

locked mass system 


and Kr, = 


2 m 2 


is tie damping ratio of the 


conpound pendulian* 

In ncjndimensional form (5*^) and (5*5) are written as 

2l ^ ^ ii |0 

p x" 4 - 2 P X* + X - p ( Q sin 0+9 cos 6 ) 




+ ( 2 C.J P ) cos t 


( 5 . 7 > 
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— n — 2 — 

and p e+ 2 C 2 P^e*+^^'“P x") sin 0 = 0 , (^. 8 ) 

respectively, where 


p = , t = at + tan”"^ (2 p) , 


(5.9) 


and the prime denotes differentiation with respect to the 
nondiniensional time t. 

Equations (5*7) and (5*8) are same as (2.4-2) and 
(2.4-3) except that the forcing function in (5*7) depends on 
the frequency and damping as well. For low damping and for 
forcing frequency near the natural frequency of the locked 
mass, i.e. p ^ 1 , the quantity under the square root sign 
on the right hand side of (5*7) is almost unity. Therefore, 
for these parameter canhinations (i.e. < <1 and p % 1 ) 

the systems, shown in Figs. 2.1 and 5*1 j have similar charac- 
teristics. Equations (5*7) and (5*8) can be solved approxi- 
mately by the method outlined in chapter 2 . 


5.3 DESIGN OF THE EXPERBIEHT 

5.3*1 Experimental Model 

The photograph of the experimental nodel is shown 
in Fig, 5*2. Fig. 5*3 shows the sectional views of the model* 
In this figin?e (a) is the primary mass, M, supported on the 
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FIG* 5*3 


S0GTIO«AI. ViaiS OF TIBS 


EypSRU^Stl'^^' MOOa 


13^ 


linear spring (b) of stiffness One end of the spring 

is fixed to the primary mass and the other end to the base 
(h). Fixed to the base is the outside guide block Cc). The 
primary mass oscillates freely within this block on strip 
ball bearings (d). The strip bearings aro free to travel 
vertically alongside the primary mass and are retained in 
the grooves as shown in the figure. On two adjacent faces 
of the primary mass, the strip bearings are pressed from 
sides by adjusting plates ig^) & (g^)* The pressure on the 
bearings can be controlled by moving these plates ty the 
adjusting screws (s). These screws aro set with an optimum 
pressure, obtained by trial and error, to ensure free move- 
ment of the priraaiy mass. Bie cempound pendulum (e) is 
hinged to the shaft (f ) which in turn ip fixed to the primary 
mass. The cempound pendulum is made of an eccentrically 
hinged disk. The effective length of the pendulum is 



(5.10) 


where D is the diameter of the disk, d is the hole diameter 


and e is the eccentricity. 


5 . 3.2 Experimental Set-up and Instrumentation 

Figures 5*^ and 5*5 show, respectively, the photo- 
graph and the block diagram of the experimental set-up with 
the associated instrumentation. The experimental model 
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shorn, in Fig» 5*2 is rigidly mounted on the table of a 
mechanical vibrator* The table provides a hamcnic dis- 
placem^t which is independent of the load acting on it. 
Both the amplitude and the frequency of excitation can be 
varied continuously and independently in the range 0 to 
0.l5 inches (total displacement) and 5 Hs to 100 Hz, res- 
pectively. An accelerometer is rigidly mounted on the top 
of block (a) shown in Fig. 5*3* 

Output from this picluip is fed to the vibration 
meter which shows the displacement of the primary mass. 
Output signal of the vibration meter is fed to a u - v 
recorder for taking permanent records and also to an osci- 
lloscope for visual display of the signal. 

The amplitude of the vibration table is set and is 
measured \/ith the help of a travelling raic rose ope. Another 
accelerometer is mounted on the vibrating table. Through a 
charge amplifier the output of the accelerometer is fed to a 
frequency counter to measure the frequency of excitation. 

The frequency is measured indirectly by measuring the time 
periods. 

The details of the instrimients used are listed 


in appendix A-5 
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5.^ OKrECTIVE AND PROGEDURL OF THE ©[PERBiSJIT 

The objective of the experiment is just to verify 
qualitatively the theoretical results of chapter 3. In ether 
words, experiments are conducted to verify the existence of 
parameters such that for a given harmonic base excitation 
the pendulum oscillates harmonically at sesae frequencies and 
sliows random oscillations at some other frequencies. 
Accordingly, experiments are conducted at certain discrete 
frequencies with a fixed base amplitude. Hie frequencies are 
chosen so as to reveal different types of motion cf the pen- 
dulum, namely, (i) steady state harmonic oscillation, 

(ii) random motion and (iii) strongly stable equilibrium 
position. Measurement of the pendulum velocity and amplitude 
is quite complicated. Since the quantitative results are 
not aimed at, the pendulum motion is not measured. Instead 
the pendulum motion is observed only visually and the response 
of the primary is recorded. 

Tlie dampings in the primary mass and in the pendulum 
are idealized as viscous type. The linear natural frequency 
of the locked mass system, and the damping ratio, 

are calculated from the records of free vibration decay. 

The spring stiffness, k^, is measured directly from static 
tests. The effective mass of the locked system, M +m, : 

including the contribution of sliding bearings, is calculated 
from the values o:^||k^ and The linear natural frequency; 
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of the pendulum, is obtained from tiie measurement of 

its time period* Since is found to be quite low, the 
number of cycles, n, to bring dovjn the amplitude of oscilla- 
tion to half of its initial value is counted and is 
calculated frcM the relation 

exp i ~ 2 n % (5*11) 

By measuring the eccentricity, e, of the compound pendulum, 
its effective length, 1, is calculated from (5.10). Inertia 
factor, is obtained from (5*6) as 

If = a,| ( i ) (5.12) 

The experimental model is given a base excitaticn 
by presetting the amplitude and vaiying tlie frequency of the 
vibration table. Responses of the primary mass and the 
pendulum are observed for vari<sus ccmbinations of the table 
amplitude and frequency. For the cases of interest, response- 
of trie primary mass is recorded. 

5.5 RESULTS AUD DISCUSS I0H3 

The measured values of the system parameters are 
as follows: 

spring stiffness =1.2 kg/cm , 

natural frequency of the locked mass system f^ = 6*4- ira, 
locked mass M + m = 728 gm, 
cempound pendulum mass m = 170 gm, 



140 

natural frequencjr of tho canpound pendiiLum f^ = 2.^ hz, 
ecc on tricity of tiio pendulum o = lA cm, 
effective length of the pendulum 1 = 1.53 cm, 
inertia factor c O.356, 

damping ratio of the locked mass system 5”^ - 0*044 
and damping ratio of the pendulum = 0.005. 

(5.13) 

With the above values the nonoimensional quantities are 
calciilated as 

R = g = 0.305 , q _ _ 0.375 j 

— y _ -1 - (5*14) 

Y = j = 0,031 and p = = ^-.4. ^ j 

^tfhere T is the time period of the base excitation in 
seconds. 

With Y = 0,132 cm ( Y = 0,031 ), over a certain 
range of excitation frequency caapletely nonperiodic respcnses 
of the primary and the pendulum are observed. Experimental 
results for I = 0,132 cm \,’itli Increasing excitation frequency,' 
are presented belov:. Theoretical responses, obtained by , 

numerical integration of (5*7) and (5.8), vith parameter ' 

values given by (5.13) and. (5.l4), are also discussed. 

Initial conditions used for numerical integrations are 

X (0) = X' (0) = e‘(o) = 0 and e (o) = 0,05 (5.15) : 

n 

and the step size of integration is taken as , unless 
mentioned otherwise. 
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Theoretical ncaidiiaensional response curve for the 
locked mass system is shown in Fig. 5-6. The circled points 
show the values obtained frcm the experiment when the pendulum 
remains stationary. 

The experimentally obtained locked mass response 
record, Xj^ (t), for the escitation time period T = 195 msec 
( p = 0,8) is shown an Fig. 5»7« In this case the pendulum 
is not disturbed and stays in its static equilibrium position. 
The response record is seen to be harmonic and the pk - pk 
value, as measured by the vibration meter, is 280 mils. This 
reading is indicated by the point (a) in Fig. 5*6. Now if 
the pendulum is distiarbed slightly from its static equili- 
brium position, its oscillations start growing until^ the 
pendulum and the primary mass reach harmonic steady states. 

The steady state respcnse record of the primary, x (t), is 
shoi'mi in Fig. 5*8 and the pk - pk value is noted as 185 mils. 
This gives the nondimen sional amplitude of the primary mass 
as 0.055* The numerical integration yields a stable harmonic 
steady state as x = O.076 cos ( t - 8.6), 9 = 0,45 cos ( ^ ‘ 

2.5). 

It may be noted that when the pendulum oscillates 
the amplitude of the primary mass is reduced to 185 mils as ■ 
ccmpared to 280 mils with the pendulum stationary. i 


f 
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By Increasing the excitation frec^uency to p = 0.82 
(T = 190 msec) similar patterns of the responses are ohserved. 
Figure 5»9 shows (t) i/hen the pendiiLum does not oscillate 
and Fig. 5.10 shows the primary response, x (t), \/hen the 
pendulum also oscillates harmonically. Vibration meter shows 
the pk - pk values for these cases as 330 mils and 285 mils, 
respectively. For the latter case the nondimensional ainplitude 
is 0.08^, Integration also gives a strongly stable steady 
state as X = 0,09 cos C t - 10. 5), 6 = 0.28 cos ( | - 0,1 ), 

The reduction in the primary mass amplitude, due to oscilla- 
ticffi of the pendulum, may again bo observed at this frequency 
as well. The amplitude of x^ (t) is indicated by the point 
(b) in Fig. 5 * 6 . , 

I 

The locked mass response record, Xj^ (t), for • 

p = 0,97 (T = 16O msec) is shown in Fig. 5*11 and the pk - pk 
value is measured as 800 mils. The locked mass response is 
again seen no be Iiamonic and the corresponding amplitude is ' 
shown by the point (c ) in Fig. 5»'3. Somewhat large discre- 
pancy between the experimental and the theoretical values of i 
the locked mass amplitudes may be attributed to the sharp 
gradient of the response curve at this frequency range. At ' 
this frequency if the pendulum is disturbed slightly frem I 
its static equilibrium position, the subsequent pendulum 
motion does not show any regular pattern even after a long 

I 

time. It oscillates irregularly for some time and then 

I 

! 
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makes few revolutions. The number and the direction of 
revolutions also do not show any pattern. Consequently the 
response of the primary, x (t), ceases to be periodic. The 
experiment is repeated at the same frequency by stopping the 
pendulum and then again letting it oscillate. After a lapse 
of sufficient time, v/hen the transients are expected to die 
down, the response x (t) is recorded. Seme such sample 
records are shown in Fig, 5»12. It can be seen that the 
response records appear different. Moreover, the records, 
though having some regularity, do not look periodic. This 
regularity in the primary motion is due to the fact that the 
primary system is directly driven by the base excitation. 

This was also observed in the theoretical responses as dis- 
cussed in section 3*5*3«^* The pendulum motion is much more 
irregular and appears to be random in nature. The randomness 
in the nature of the pendulum motion is attributed to tlie 
inherent external disturbances and the fact that the para- 
meters such as the excitation frequency etc, do not remain 
exactly at the same value throughout the experiment. Even 
if there exists a periodic steady state, it is so highly 
sensitive (for a particular combination of parameters) to 
those random disturbances and cimngos in parameter values 
that the pendulum can not attain this state and hence osci- 
llates randomly. At this excitation frequency the numerical 
integration, with the initial conditions given by C5-15), 
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also does nof converge for decreasing step sizes h = |~ , 
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At the excitation frequency p = 1.01 (T = 155 
msec) similar random motion of the pendulum is observed# 

An increase in the frequency p to 1*12 (T = ilfO msec) makes 
the static equilibrium position of the pendulum stable. Iho 
disturbances to the pendulum die out and the system oscilla- 
tes like a locked mass with a pk - pk value as 275 mils. 

The response is found to be harmonic and the corresponding 
aJi^jlitude is shown by the point (d) in Fig, 5*6, 



CHAPTER 6 


AUTOPARAMETRIC VIBRATION ABSORBERS 
6,1 D^TRODUCTICf^ 

Haxton and Barr [ 20 ] studied the possibility of 
using an autoparametric vibration absorber where the parame- 
trically excited secondary system neutralises the vibrations 
of tho primary system. In this chapter some observations 
are made cxi the range and scope of the absorber acticxi of 
tho pcaidulum on the primary mass for the system studied in 
previous chapters. The conventional auxiliary mass damper 
is also studied briefly with an additional degree of froedem 
in rotation provided to tho ainciliary system. The natural 
frequency of this rotational motion is adjusted so as to cause 
internal resonance, 

6.2 AUTOPARiLMETRIC VIBRATION ABSORBER (TWO DEGP:® - OF - 
FREEDOM SISTM) 

In the system considered in chapter 2 (Fig. 2.1 ) 
the pendulum is excited parametrically. With a proper seleo« 
tion of parameters tho pendulum acts as an absorber for the 
primary mass as is evident fron Fig- 2,6, She parameters 
which affent tiic ansorber characteristics are tho mass ratio 
R, the damping ratio of the pendulum ^2 frequency 

ratio q. The effect of changing q is shown in Fig, 6,1. 






All the rc cults presented in Figs. 6,1 to 6,^ are 
obtained bv the second order approximation in harmonic balance. 
These results match well -with the integrated results except 
for small differences in the overhang regions. 

In Fig. 6.1 tlie amplitude of the primary is expre- 
ssed as the ratio A/P i.e. the dynamic magnification 
factor. Curve 1-1 is obtained with = 2 » ^2 ~ ^ 

i.e, the linear natural frequency of the pendiilum ( “2^ 
taken as half of that of the locked mass. Curve 2-2 is 
with q^ = I - , q2 = 0 i.e. “2 natural 

frequency of the primary mass. The location and the magni- 
tude of the peak are seen to change with change in q. 

Curve 3-3 shows the response of the primary astern, witti 
the same excitaticn, in absence of the auxiliary system. 

Figure 6.2 shows the dynamic magnification factor 
for different values of The response in the overhang 

region are indicated by dotted lines. It can be seen that 
(with the combination of parameters as noted in the figure) 
the optimum absorber action is obtained with ^2 — 0,10, 

It was shown in chapter 2 that vltb. ?2 = 

(R = 0.2, = 0.02, q^ = 0,5), over a certain frequency 

range, tiio responses become nonperiodic with P = 0,05. In 
this case if the damping in the pendulum, ^2, is increased 

to 0,15 the responses are again harmonic and are shown in 
Fig. 6.3* 



FIG, 6*3 AMPLITUDE (F THE PBHDDLOM AHD DBIAMIC 
KAQIIPICAIIOE FACTCa TOBSUS FRBQinKCI ; 

p a 0*05, % » 0.15 
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The locked mass rcsponsG is also included in Figs* 
6*2 and 6 . 3 . Comparison of those figures ( = *^»15) 

reveals that the dynamic magnification factor is not indepen- 
dent of P (i.o. the response of the primaiy is not linearly 
dependent on P) and the range over which the pendulum affects 
the primary response is also different for different values of 
P, Dependence of A ai P not being linear 'was also observed 
in Fig, 2*k-, 

It can be seen from Fig. 2.^ that for P < O.OO18 
the internal rescnance does not take place and the systcmi 
responds like a locked mass. For 0,00l8 < P < 0,005 the 
primary respense is almost independent of P, But for these 
small values of P the primary unstable region is also quite 
narrow (0.875 < P < 0,96 for P = 0,005 and = 0.05 vide 
Fig, 3 . 3 ) and only near the resonant frequency of the locked 
mass system (i.e, near p = 0,913) the pendulum reduces the 
amplitude of the primary. As P is increased beyond 0,005 
only the initial change in A is somewhat linear (Fig, 2,4), 
Moreover, variation of A with P is obviously dependent on 
p which is clear from Fig, 6,4, In this figure changes of 
A with P are shown at two different frequencies, Iho 
chain dotted lines show the locked mass response amplitudes 
over the range of P vhere the static equilibrium position 
of the pendulum is stable. The dotted lines show the 
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amplitudes in the OTorhang regions and tho firm lines show tho 

amplitudes v/hon e also' oscillates hamiionically# Considering 

also the fact that the responses boccxnG nonperiodic ofver 
a certain frequency range for values of P beyond a 
certain critical value (which obviously depends on 

is clear that for different ranges of P different values of 

^2 will provide the optimum absorber action. 

6.3 AUTOPARAMSTEie VIBMTICN ABSOEBER 

(TEJBBE BSGEEB - OP - PEEEDa4 SISTM) 

6.3*1 Eqmtions of Motion 

Figure 6.5 shows a threo-dogree-of-fruedom system 
■sdiich is similar to that shown in Fig. 2.1 except that the 
pendulum (under gravity field) is now considered to be 
elastic. The stiffness of the elastic pendulum is repre- 
sented by k^, c^ being the damping coefficient. The length 
of the pendiilum under static equilibrium is 1 and z repre- 
sents the relative displacement betwuon M and 

The equations of motion are 

(M + m) X + k^ X + c^ X + m z cose - 2 m z* q sine - 

m (1 + r) ( e sine + e'^ cose ) = F cos ut, (6#l) 

m z + k-, z + C-, z’ + m X cose - m (l+z)6^+mg (1- 
3 3 

cose ) = 0 , (6.2) 

2 - 

m(l+z)9 +C 2 6’ + 2ni(l4-z)ze +(l+z)(mg- 

m X ) sin 6=0 


(6.3) 
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In non dimensional form these equation are written as 


2 r-. 


(1 + R) x" + 2 p x’ + X + R p [ z” cos S - 

2 z‘ e' sine - (1 + z) ( S" sin e + 9* cos e ) ] 

= P cos t , (6.4-) 


p'^ z'* + 6^ z + 2 |i p z* + p^ x” cos e - (1 + z) p"" + 

(1 - cos 0 ) = 0 , 

2 

C 1 + z) p^ 9 " -f- 2 Cp i p 0 * +2 (1 + z) p^ z* e ’ + 


2 


(1 + z) ( q^ - p^ x”) 


sine = 0 


(6.5) 


( 6 . 6 ) 


respectively, Ihe various ncaidlmensional parameters are 
defined as follows: 


X=^jZ=^,t=(ut, fl 


1 » “ 2 "J 1 ’ 


Q 


if - 


_ /is 


V m 


> P = 


(I) 


6 = 


n. 


n. 




q = 




? _ 1 — 

1 ~ 2 M ’ 


'2 = 




P = 


1 • 


(6.7) 


With high values of 6 and p, the pendulum behaves 
like a rigid one ( z = 0) and (6.4) and (6.6) are same as 
( 2 . 7 ) and ( 2 . 8 ). If the rotational motion of the mass m 
is locked ( e(t)=0¥t) the system bdiavcs like an 

ordinary auxiliary mass damper and the equations of motion 
bee one 


p^ (1 + R) x” + 2 p X* + X + R p 


2 .7 II 


P cos t 


( 6 . 8 ) 


I 
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2 •*- 2 “• 2 — 
and p z”+6 z+2 [ipz*+p x” = 0 


(6*9) 


mie responses for the primary and the ainciliary masses are 
obtained from (6*8) and (6*9) as 


s = A cos ( t + 9^ ) 
and z = Z cos ( t + <5^ ) j 


•fehere A » Js^ + s^ » 9 ^ = tan"’^ ^ ^ 


1 

= > 9' ^ ^ - 1; ^ 


3 


( i = 1, 2, 3, being given by 

( 6 ^ - p^) + (2 p) s^ 


2 ^ o2 

s^ + Sg. 




( ^^ - p^ ) " ( 2 p) s^ 


2 2 
^6 


’3 


and S' 


'if 


P p 

“2 2 
^5 


P Sg p 

2 2 » 
S5 + Sg 


( 6 . 10 ) 


(6.11) 


vith Sr 


and s 


6 ~ 


(1 - p^) ( 6^ - p^) - (2 p) (2 p p) - R 6^ p^ 

(1 -p") 2 pp + (6^-p^) 2 c^p-p^R (2|i.p)- 

( 6 . 12 ) 
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N 


For = 0 the optimm values of 6 and ^ are 


6 = 


1 + H 


and 


(6.13) 


V 8 (1 + 


E)' 


6.3.2 EESULTS AND DISCUSSIONS 


Equations (6.^) to (6.6) are integrated ntaaerically 
witii E = 0.2, P = 0.02, = 0.02 and Initial 

cooditloKis 6 ( 0 )= e *( 0)=0 yield the results for the cen- 
v'entional auxiliary-itiass-dangicr, which can he obtained 
directly fran (6.10) to (6.12). For this case the qptiatuia 
values of g and {i (with E = 0.2, 5^=0) are obtained from 

(6.13) as 0,833 and 0,208, respectively. With these values 
of 6 and p the (dynamic magnification factor of the primary 
mass is shown in Fig. 6,6 (b) ty the chain dotted line. 


It is well known that for the values of 6 and p. 
different from the optimun values, the locations and the 
magnitudes of maxima and minima in the primary response curve 
change. For example, the dotted lines in Fig, 6,6 (b) shows 
the dynamic magnification factor for the primary mass with 
6=1 and p- = 0,1 (e(t)=0¥t). With this change 
the amplitude decreases near p = 1 but the peak at p = 0,8 
becomes quite large. It is obvious that if the auxiliary 



FIG* 6#6 (a) JiMPLITODE OF PM2XJLCM 7ER8US EXCITATIC3N 

FSEQUMCI s p - 0.02 
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mass is given an additional degree of freodem in rotation j 
its effect on the primary response will he different i To 
study this effect near p = 0*8, the. elastic pendulum is 
tuned so as to cause internal resonance near this excitation 
frequency • 

To this ^d q is set as O.if. The integrated 
results are obtained with 0 (0) = 0,1 (all other initial 
conditions are zero). These responses are found to be har- 
monic, The firm line in Fig, 6.6 (b) shows the amplitude of 
the primary whai the pendulum also oscillates harmonically 
(i,e, the static equilibriim pcsitiai of the pendulian is 
unstable). These aii5)litudes of the p^dulum moticn are shown 
in Fig, 6,6 (a). Outside the primaiy unstable region 
(p < 0,735 and p J> 0 , 89 ) the static equilibrium position 

of die pendulum is stable and the system respends like an 
auxiliary mass damper. Because of the high damping in the 
pendulum, the overhang regions are almost nonexistent. The 
pendulum notion is seen to suppress tho peak at p = 0,8, 

The amplitudes of the relative displacement of 

z 

the auxiliary mass (expressed as ^ ) are shown in Fig, 6,7, 
In this figure the representation of each curve is same as 
in Fig. 6,6 (b). It may bo noted that detuning of 6 and 
p , from their optimum values, results in a high peak near 
P = 0.8 which again is suppressed by providing the additional 
degree of freedom ( 8 - coordinate) to the auxiliaiy mass. 







CHAPTER 7 


CONCLUSIONS 

1 

This work concentrates mainly on a two degree - 
of - freedom autoparametric syst^i under a harmonic exci- 
tation. The major objective is to detemine the regions 
and nature of the harmonic and nonharmonic responses. The 
method of harmonic balance is used to obtain different 
order approximate results and these results are verified by 
numerical integrations. The major conclusions of the work 
are listed below. 

1. The method of harmonic balance is sufficient to predict 
the existence and stability of tiie harmonic or near 
harmonic steady states. The amplitude and tlie stability 
conditions are obtained in closed form only for the 
first order approximations. 

2. dily for sufficiently small excitations, the first 
order approximation (which essentially solves the 
approximate equations of motion with quadratic non- 
linearities) yields correct results. 

3. The higher order approximations re-'^'eal that (i) the 
primary amplitude is not independent of the excita- 
tion amplitude and Cii) the stability behaviour 
predicted by the first order approximation may be 
incorrect. 
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r’' 

For tile system having a spring controlled pendulum, 
v;ith certain critical values of the amplitude and 
frequency of excitation, tlie stable harmonic motion 
bifurcates to one unstable harmonic motion and another 
stable amplitude nodulated periodic state. Ilie bifur- 
cation is associated with a jump in the amplitude 
values# 

5# For the system having a gravity controlled pi^dulum, 

, with certain critical values of the forcing ampli- 
tude, foLloifing two types of brandling (depending on 
tte frequency) of the stable harmonic moticn are 
possible: 

(i) With a moderate excitation the stable har- 
monic motion branches cut to one (or more) unstable 
harmonic state (s) and one (or more) amplitude nodulated 
periodic motion(s). In this branching the amplitudes 
do not increase sharply. 2he amplitude modulated 
suates are weahly stable and the modulation periods 
can vary over a wide range. 

(ii) Ihe stable harmonic motion branches out to 
one (or more) unstable hai**ionic state (s) and another 

cliaotic state. The transition from the stable harmcnie 
motion to chaotic motion is associated v/ith sharp 
increase in the amplitudes. 
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6. For the conhinations oi' forcing amplitude and frequency, 
resulting in chaotic states, the numerical integrations 
do not converge and the responses are completely non- 
periodic though remaining hounded. 'In such cases a 
meaningful description of the integrated responses is 
'possible throu^ the statistical analysis. 

7. Randan (or nonreproducible) responses are obtained 
experimentally xcLth certain canbinations of the forcing 
frequency and amplitude for which the numerical inte- 

. grations also do not converge. 

8. The pendulum in the two degree- of-freedon syston exhi- 
bits limited absorber action. The major disadvantage 
of this absorber is that idle effective frequcaicy range, 
the level of the neutralizing effect and the optimum 
value of the absorber damping depend upon the excita- 
tion amplitude. 

9* With a nonop thnal auxiliary mass damper, the peaks in 
the responses may be suppressed by providing the 
secondary witi an additional degree of freedom in rota- 
tion (so that it behaves like an elastic pendulum). 
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APPENDIX A«1 


DEFINITIONS USED IN STATISTICAL ANALYSIS PRESENTED 

IN CliAPTBR ? 

Various definitions and relations used in tho 
statistical analysis are taken from reference [48] • 

A single time history representing a random phono- 
menen. is called a sample function (or a sample record when 
observed over a finite interval). The collection of all 
possible sample functions which the random phenomenon might 
have produced is called a random process or a stochastic 
process. 

For a given random process {u (t)) , where the 

symbol { } is used to denote an ensemble of sample functions, 
the mean value (t^) and the autocorrelation function 

\ ti + X ) can be found by taking the ensemble average 

at time t^, Wlien ^u ^ vary 

with t^, the random process {u (t) ) is said to be weakly 
stationary. If all the possible moments and joint moments are 
time invariant, the random process is said to be strongly 
stati alary. 

If for a given randan process the time - averaged 
mean value and autocorrelation function are equal to the 
corresponding ensemble averaged values, the process is said 



to bo a weakly orgodic rando,.; process. For n strongly 
ergodic process all tiio onsoinblc averaged statistical proper- 
ties are doduciblo from corresponding time averages, k samplo 
record obtained from an orgodic process will bo stationary. 


APPENDIX A- 2 


BASIC DESCRIPTIVE PROPERTIES OF A RANDOt^ DATA 

For an orgodic randan process Cu (t)} various 
parameters arc defined as followsJ 

(i) Moan value is given oy 

1 ^ 

T / u dt (A-2.1) 

T-. CO - ° 

♦ 

p 

(ii) Variance is given by 

T 2 

^ ~ Lt» J' (. u (t ) •“ ^u ^ CA’-’2#2) 

rp O 

X ->1 oo 

(iii) Autocorrelation function R^ ( t ) is given by 

i T 

R^ ( ) = Lt. ^ j u (t) u (t + ) dt , 

rp O 

(A-2.3) 

where 'c is the tine displacement. 

R^ ( T ) is a real valued even function with a 
maximum at r = 0. A transformed deterministic data 
with zero mean has R^ ( '^ ) which persists over all 
T , as opposed to the randan data for which R^ C ) 
approaches zero. Thus a study of an autocorrelogram 
can reveal the deterministic trends in apparently 
random data. 
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(iv) Po'uer spectral density (psd) function establishes the 
frequency composition of the random data. One sided 
PSD function can be defined as the general frequency 
canposition of tlio data in terms of the spectral den- 
sity 01 its mean square value. By this definition PSD 
is evaluated directly from the data points. Alternate 
dsfinition of PS!D function is in tcpiiis of FoupiGP 
tPcinsforrA of tlio nut oc or relation f unction ♦ Iliis gives 
a two sided PSD function as 

CO 

^u \ ^ ^ (-2i7i f r ) a , 

— oo 

(A-2A) 

whore f is the frequency at which is being 

evaluated and - oo < f < ■» , 

For f varying over positive range only, the one 
sided PSD function is 

oo 

(f) = if / \ ) COS (2ti f T ) dT *, 

o 

o ■%. I < ^ , (A-P-5'> 



APPENDIX A -3 


RUN TEST TO CHECK THE HIPOTHESIS OF STATICDTARin OF 

A SAMPLE RECORD 

The stationarity of a randcai process is tested try 

investigating a single record u (t) as follows: 

(i) The same record is divided into J equal time inter- 
vals where the' data in each interval may be considered 
independent. 

Cii) Variance, is computed for each interval* Omitting 

the subscript u , these values are aligned in time 

^ 2 2 2 2 
sequence as cr ^ i ^2 » c^3 » ••• j oj • 

(iii) Finally this sequence is tested for the presence of 
underlying trends or variations other than those due 
to expected sampling variatiens. This is done by run 
test. Run test uses a nonparametric approach which 
does not require a knowledge of the sampling distri- 
butions of data parameters. The rim test is perfoimed 
as follows: 

Let it be hypothesized that each variance in the 
sequence ( cf-jj independent 

sample value of a random variaele. If this hypothesis 
is true, the variations in the sequence of sample 
values will be random and display no trend. Hence the 
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niJimber of runs in the sequence relative to, say, median value 
will bo as expected for a sequence of independent randcin 
observations of the random variable, which is presented in 
the Table in appendix A-4, If the ninnber of runs is signi- 
ficantly different from the expected number given in the 
Table in appendix A-4-, the hypothesis of staticaiarity would 
be rejected. 



Table for Percentage 


Values of r„. such that Proh 
n, a 

1 = Number of observations 


4 i 5 

^ 5 0.99 § 0,975 5 0.95 


10 

5 

6 

6 

i5 

9 

10 

11 

20 

13 

14- 

i5 

25 

17 

18 

• 19 

30 

21 

22 

24- 

35 

25 

27 

28 

4-0 

30 

31 

33 

^5 

34- 

36 

37 



50 


38 


4-0 


S A-4- 


its of Run Distribution 


( r > r ) = a , 'where n = 

ii. 5 CC 



26 27 28 

32 33 34- 

37 39 4-0 

4-3 44- 4-6 

48 50 51 

54- 55 57 

59 


61 


63 


rol'^ 



APPENDIX 


DESCRIPTIONS OF THE INSTRUMENTS USED IN THE EXPERIMMT 

1. Vibration Fatigue Testing Machii 
Model 150-VP-T , (Vertical) 

All American Tool & Manufacturing Caapany 

2 m Vibration Picloip 
Type I 5 &O - P52 
General Iteidio Coni^ny 

3» Vibration Meter 
Type 1553-A 
General Radio Conpany 

4. Visic order Oscillograph 
Model 9O6T - 179OFOO 
Honey Well 

5« Columbia Model No. 502-1 Ac c el er erne ter 

6. Electrostatic Charge Amplifier 
Model 566 M3 Multi- range 
Kistlor Instrument Gorporation 

7. Digital Universal Counter and Timer 
Model - 613 



Yamuna 



